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MSP SPECIFICATIONS 


OBJECTIVE: Return IKg martian soil sample from the surface of mars 
to a mothership in a 600km mars orbit. 


STRUCTURE: 

overall height : 1.2 meters (touchdown configuration) 

1.5 meters (stowed/descent configuration) 
overall width : 0.66 meters (touchdown configuration) 

1 .5 meters (stowed/descent configuration) 
maximum gross vehicle mass: 110.25 Kg (after mothership 

separation) 

mass at return: 2.00Kg (includes IKg sample) 
planform: 3-legged layout with hexagonal bus 
materials: Beryllium(launch rings, sample retrieval system, 
heat shield backing, landing gear main struts) 
Magnesium (shock absorbers, landing gear secondary 
struts, bus frame) 

Aluminum (bus platform, vehicle skin) 


PROPULSION 

deorbit/descent/ascent: 1 single restartable liquid propellant 
motor, pressure fed, using 50%UDMH, 50%Hydrazine 
as a fuel and nitrogen tetroxide as oxidizer 
apogee kick(rendezvous): solid propellant motor 
maneuvering: (8) thrusters using UDMH as propellant 


MISSION PROFILE 


This mission begins with detachment from the mother ship, 
the Mars Orbiting Vehical (MOV) of the designed spacecraft, the Mars 
Surface Probe (MSP). After separation the MSP determines its 
position and the relative position of the landing sight and initiates a 
short burn which will cause the MSP to be at our optimum position 
for deorbit within seven days. During this time the ship is in a 
dormant mode with occasional systems checks made by a monitoring 
microcomputer. Adjustments to the precession orbit are 
periodically made. When the MSP is in the correct location to 
descend with a minimum fuel expendature and no plane change, t e 
main motor fires, decelerating the ship. 


The MSP turns around to face its heat shield toward the 
direction of flight. When the atmosphere is encountered at about 
200-1 50km altitude the ship releases a kind of "drogue chute for 
passive stability during descent. 

2-3 

The MSP sheds velocity during reentry, protected by its 
heatshield and a multilayer insulation. At 1km above the actual 
ground surface a parachute is deployed to slow the craft even 
further. The heat shield is jettisoned with explosive bolts. After 
slowing to 42 m/s and drifting to within 100 meters of the surface 
the chute is released and the main motor fires to bring the craft to a 
standstill 5 meters above the surface. 

3 4 After stopping 5 meters above the surface the main motor 
shuts down and the MSP falls freely to avoid contamination of soi 
samples. Once safely on the ground a shovel type mechanism deploys 
and begins collecting sample soil. 

4 5 Once the sample has been collected and stored the ship waits 
for half a martian day until the MOV orbit will be nearly overhead. 
The MSP ascent stage separates from the hexagonal descent stage 
leaving spent fuel tanks, the sample collector, the landing gear, and 

other used parts. 

The ascent stage climbs out of the atmosphere using the mam 
liquid motor for the last time. 



The ascent stage burns out. The computer points the rendevous 
stage in the right direction, spins it up and releases it. There is no 
active guidance on the rendevous stage. 

7 - 8 

The rendevous stage makes its solid motor burn and activates 
its beacon for the MOV to find it. Its orbit differs slightly from the 
MOV so the two will eventually pass each other. 

8 - 9 


The MOV catches up to the rendevous stage and brings it inside 
the docking bay which originally held the full MSP and then secures 
it for the return trip. 



7: Solid apogee motor burn and orbital insertion 




Descent through atmosphere, deployment 
of guide/drogue chute 



MSP Mission Scenario 
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3; Chute release, main engine retro firing 
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5 : Ascent stage separation 


4: Sample acquisition 
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Mass: 


Final (on orbit) : 

Dirt: 

Sample Container: 

Apogee Engine Case: 

Beacon: 

Batteries/solar cells: 
Structure: 

Total : 

Before Separation: 

Apogee Stage: 

Apogee Fuel: 

Apogee Stage Support: 

Fuel and Helium Tanks: 
Valves and Regulators: 
Aeroshell : 

Guidance and Control : 
Batteries: 

Liquid Engine and Plumbing: 
Engine Structural Support: 
Thermal Insulation: 

Attitude Control: 

Total : 


Ascent Stage Before Liftoff: 

Ascent Stage Dry Mass: 
Ascent Fuel (5% ullage) : 
Total : 


Mass on Mars: 

Ascent Stage: 

Descent Tanks and Plumbing: 
Landing Gear: 

Payload Scoop: 

Batteries : 

Laser Range Finder: 

Landing Computer: 

Attitude Control: 

Thermal Insulation: 
Structure: 

Total : 


Mass at Re-entry: 

Landing Mass: 

Retro Fuel : 

Parachute : 

Parachute Canister: 

Heat Shield and Supports: 
Drogue Basket: 

Total: 

Initial (On-Orbit) Mass: 

Re-entry mass: 

Descent Fuel : 

Orbital Maneuvering Fuel: 
Total : 

Table ( 

MSP Mass Estimates 


1.00 (kg) 
0.10 (kg) 
0.25 (kg) 
0.30 (kg) 
0.15 (kg) 
0.20 (kg) 

2.00 (kg) 


2.00 (kg) 
0.73 ( kg ) 

0.90 ( kg ) 

3.10 (kg) 
0.50 (kg) 
1.80 (kg) 
2.00 (kg) 
0.10 (kg) 
2.00 (kg) 
0.90 ( kg ) 

0.50 (kg) 
0.50 (kg) 
15.03 (kg) 


15.03 (kg) 
44.17 (kg) 
59.20 (kg) 


59.20 (kg) 
1.50 (kg) 
1.50 (kg) 
15.00 (kg) 
0.80 (kg) 
1.00 (kg) 
1.00 (kg) 
0.50 (kg) 
1.50 (kg) 
7.55 (kg) 
89.55 (kg) 


89.55 (kg) 

3.00 (kg) 

2.00 (kg) 
0.30 (kg) 

5.00 (kg) 

2.00 ( kg ) 
101.85 (kg) 


101.85 (kg) 
7.57 (kg) 
0.83 (kg) 
110.25 (kg) 
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Section A: 

Structural Design and Component Placement 


Eric Deyerl 


1 . Goal 


To provide the most compact (space and weight efficient) structure to unite, 
protect, and allow effective operation of all spacecraft components. 


2. Development 


2.1 Research 


Much research was accomplished, mainly through the perusal of large amounts 
of NASA technical reports (see References). Here, the desianer focused his 
attention on reports that dealt with landing gear and spaceframe studies. Most 
specifically, there were two NASA projects that provided the most insight in 
these areas; those were the Surveyor and Apollo programs. (Figures 1 ana 
2).Although both of these missions were lunar in objective, they had similar 
parameters in their missions, and so were valuable in the information they 
provided. 



Figure 1 : Surveyor spacecraft 



Figure 2: Apollo lunar lander 



Figure 3: Generic three- 
legged lander 
planform 


The surveyor three-legged layout was inspirational for the MSP designer. It is 
simplistic and effective, and nad been used successfully in the Viking mission to 
Mars. (Refer to Figures 3 and 4) With three leas, one can be sure that all 
footpads will be in contact with the Martian soil upon landing. The three-legged 
design also allowed nicely for a structural-rigid hexagonal bus structure. Finally, it 
is also a stable platform whose shape makes it difficult for surface winds to 
catch a broad side of the craft and topple it. 


Figure 4: Viking spacecraft 


PASE IS 

QUALITY 





The surveyor project also gave valuable insight into landing gear and shock 
absorber design, in conjunction with the placement of crushable honeycomb 
material. (Figure 5) 
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Figure 6: Apollo lunar lander 
landing gear leg 


The Apollo Lunar Module's landing gear was also studied in an effort to 
incorporate a similar retracting system on the MSP. This was later found un- 
necessary. A rendition of one of the "Eagle’s" gear ,s found in Figure 6 . 



2.2 Preliminary Layout Studies 

Prior to much of the research, preliminary sketches were made to seek out 
viable spacecraft layouts. Two early plans and some notes on each are given in 

Figures 7 and 8. 
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Figure 7: First preliminary MSP layout sketch and comments 
Figure 8: Second preliminary MSP layout sketch and comments 
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Later, after more propulsion studies and relative sizing of components had been 
undertaken, it was possible to produce more accurate and intermediate layouts, 
as shown in Figure 9. 



It was at this point that the designer realized that the critical components, as far 
as space was concerned, were tne fuel and oxidizer tanks, (which is usually the 
case with light aerospace vehicles). Thus, a fuel tank study was initiated, which 
(as shown in Figure 1 0) resulted in a possible toroidal tank configuration. This 
idea was later discarded due to complexity of purging, pressure containment, 
and the like. 


CYLINDERS 



TOROIDS 


Volume- r h - 17,200 cc 


r 

h 

10 

55 

15 

24 

20 

25 

25 

9 



r 



Volume - [ ((r - r )/ ) ][ ((r + r )/ )] 
(Approximation based on r - r ) 


r 

r 

8 

26 

10 

25.5 

12 

25.25 

15 

24 


Figure 10: Fuel tank study 






2.3 Final MSP Configuration 

After many hours of conference, deliberation, and improvisation, the layout of 
the MSP was finalized to that shown in Figures 1 1 and 1 2. 

It is a space efficient design, and yet tries to remain as simple as possible. The 
design is roughly axisymmetric, with the ascent vehicle able to clear all bus 
components. 


3 Component Design 

Following are the steps taken to design the critical components of the MSP, 
which include the landing aear, shock absorbers, bus frame and plattorm, and 
ascent and apogee stage launch rings. 

Care was taken to use the lightest and (naturally) strongest materials possible, 
with thought given to their manufacture and interaction with the environments 
the MSP would be subjected to. It also was always desired to use a safety 
factor of at least 1 .5. As will be shown, because of the generally small 
dimensions of the craft, this would never be a problem. 

3.1 Landing Gear 

It was chosen to use the landing gear configuration shown in Figures 1 3 and 1 4, 
with a main strut containing a shock absorber to cushion the vehicle at landing, 
and secondary struts for any torsional forces encountered by the gear. The 
following are the calculations made in the design of these components. 

3.1.1 Loads on Gear 

The mission scenario calls for the MSP to impact the Martian surface at 5.5 
m/s. The force on one gear (considered in the extreme case of a one-legged 
landing) is given by 

Figure 15: One-legged landing 

F = ma 

where m = 73.9 kg (at landing) 

a = (v2/2s) 

where v = 5.5 m/s 

s = 0.20 m (vertical gear travel) 

Thus, the force on one leg (taken to be applied through the axis of the shock 
absorber in the case of a lopsided one-legged worst-case landing, depicted in 
Figure 1 5) is 

Landing force on one leg = F = 5589 N 





Figure ll.uMSP Exterior Dimensions 

(Descent Configuration) 






Figure 1 1 : MSP Cutaway 

(Descent Position) 
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Figure 1 2.2: MSP Planform Exterior 








Figure 1 2.1: MSP Planform Dimensions 




Figure 1 3: Landing Gear Geometry 



1 4 cm 


Undeflected Gear 

(Before Landing) 







3.1.2 Shock Absorber 

The shock absorber chosen was a fluid-filled type, in which the travel of the 
piston forces fluid through a small orifice in the piston, damping the spacecraft s 
motion. A cutaway of the MSP shock is given in Figure 16. 


Figure 1 6: Shock Absorber Cutaway 


8.4 cm 



From the derivation given in class, which uses Bernoulli's principle to equate mass 
flow generated by the motion of the piston of area A to the resulting mass flow 
through the orifice of area A*, the following relation results: 

A/A* = (2F/pA) 1/2 (1/v) 


w h ere F = 5589 N (landing force) 

p = 917 kg/m 3 (assumed from SAE 30W oil) 
v = 5.5 m/s (landing velocity) 


which results in 

A* = (1.57)A 3 / 2 



Then, for the sake of a compact design, choosing a piston radius of 1 cm gives 


~~ Piston area = A 3.14 x 10'^ m 2 ~ 

Orifice area ° A* ■= 8.74 x 10*^ m 2 

From the geometry of the landing gear, it can be seen that to effect a 20 cm 
vertical deflection of the footpad, it is necessary to have 

Piston stroke travel = 8.4 cm 

To design the dimensions of the shock cylinder for strength, it was necessary to 
find the maximum pressure developed in the cylinder. Here, it was assumed mat 
immediately upon impact, the maximum force of landing was transmitted to the 
shock absorber (this is obviously a more-stringent-than-reality assumption). So 
here, the pressure is given by 

P = F/A 

where again F = 5589 N 

A - 3.14 x 10‘ 4 m 2 

Maximum pressure in shock cylinder = 43.76 atm 
(44.36 x 10 5 Pa) 

Finally, to find the maximum stresses occurring in the cylinder, we use the 
relations from Ugural, 

a r = (rj 2 )(P)(l - (r 0 2 /rj 2 V(r 0 2 - r i 2 ) 

o e = (n 2 )(P)(l + (r o 2 /r; 2 ))/(r 0 2 • rj 2 ) 
a z - F/(p )(r Q 2 - rj 2 ) 

where here, rj = 0.01 m 

r Q = 0.01 4 m 

p » 44.36 x 10 5 Pa 

This results in the following 

Maximum stresses in the shock cylinder 


<r r = 4.44 MPa 
a e = 24.59 MPa 


As for the material chosen, at first Lockalloy, the Be-Al metal looked promising, 
but proved too heavy for original gear weight estimations, so the following metal 
was chosen (Bruhn): 

' Material for Landing Gear: AZ61 A Magnesium ™ 

Ultimate Strength: 96.53 MPa 

Density: 1791 kg/m® 

Obviously, our safety factors are well-satisfied, and as long as the Magnesium 
can be machined to 4 mm, there should be no problem. In fact, the weight 
savings are tremendous: 

Mass of shock absorber: 

Fluid = (stroke)(area)(density) = (0.084)(0.000126)(917) = 97.1 g 
Cylinder = (length)(area)(density) = (0. 1 0)(TT )(0.01 2 - 0.01 )(21 00) 

= 29.0 g 
TOTAL = 0.189 kg 

In addition, the main strut must be checked for Euler buckling, where the critical load 
applied to the strut is given by 

P cr = 7t 2 EI/L 2 

After performing this analysis once, it was found necessary to use Beryllium, where 

Material tor Landing Gear Main Struts: Lockalloy (Be 
Young's Modulus: 2.9 x 10^ Pa 
Yield Strength: 3.24 x 10® Pa 

Density: 1826 kq/m3 

E = 2.9xl0 ]1 Pa 

| = nr 4 /4 = (7t)(0.01 ) 4 /4 = 7.85 x 1 0 * 9 m 4 
L = 0.40 m 


And thus 

Critical buckling load for main strut = 1 40,400 N 

Which is well within the factor of 3 for safety used for most buckling and instability cases, 
considering that the applied load will be about 5900 N. 


Mass of gear struts: 

Braces (2) = (2)(length)(area)(density) = (2)(0.335)(lf )(0.005 2 )(1 791 ) 

- 94.2 g 

Main Strut = (0.422)(ir)(0.01 2)(1 826) = 229 g 
TOTAL/leg = 0.323 kg 

Total mass of landing gear = (0,0189 + 0.154 + 0.09421(3 ) + footpads = 1 50_kg. 







Figure 1 8: Analytical Modelling of Bus Shelf 




3.2.1 Platform Design 

The bus platform .encounters tK^ing^af its 
inertial force of the a ? ce ^ v ^ ,C I e 5 ^ 0 fHs hexagonal shape, the platform 

r: modelled^ ocular plate under a distributed loading, as shown m F.gure 
18. From Urugal, we have 


°rmax 


= (3/4)(P 0 )(r 2 )/(t 2 ) 


r\ 

where P 0 = load/area = (m asc )(a)/(7t)(r ) 

= (45kg)(75.6m/s 2 )/(x)(0.32m 2 ) 

= 33.2 KPa 

t = thickness of plate = 0.025 m 


Thus Chnax = (0.75)133, 200)/l TC ) = 7.64 KPa 


From Bruhn, we are able to find a honeycomb sandwich material as shown 
Figure 19: 

Material tor plat 

Flexural rigidity: approx 4.5 MPa 
Density: 76.4 kg/m^ 


So, for our needs, the platform will have a mass of 



And obviously, well within the range of safety. 
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Figure 19: Platform material construction 




Figure 20 : Analytical Modelling of Bus Frame I-Beams 



0.5 cm 




From the classical beam theory analysis, the maximum stress occurs from bending 
at midpoint in the beam, and is given by 

Gb = M(h/2)/l 

Where M *> moment applied *= (landing force)(half of beam length) 

« (5589N)(0.20m) = 1 1 1 8 Nm 
h ■ height of beam = 0.20 m 

And where I for the beam is calculated by taking the moments of area each of 
the portions in Figure 21 separately, as is classically done. So, using the 
notations in Figure 21 , we have b 


Figure 21: I-Beam Cross-Section Dimensions 



Portion 


Area 


Vrm — ) 


Own i 


(l)(b) |0.5)(a+t) (l)(b)(l/4](c + t) 2 {b)(l 3 )/l 2 

(l)(o) 0 0 (l)(a 3 )/l 2 


same as portion 1 


Inserting all the relevant dimensions gives 
I section = ^ .61 x 10 ^ m 
Which ultimately leads to 

Maximum stress in beams = 1 4.7 MPa . — 


Which is well below the same ultimate stress level of the Magnesium. As for the 
masses of the beams, we can use that 

Mass = (cross-sectional area)(length)(density) 

I [\ 2 H0 C . 005 ) (0^05 )+(0 > 1 9)(0.005)](0.40)( 1 791 ) 

= 1 .04 kg/ beam 



3.3 Launch Rings 

tffZ? K, s- teep^e^ra ft together 

the ^ major vehicle forces of chute deployment and landtng. The two 
rings are depicted in Figure 22. 


Figure 22: Launch Rings 



4.75 cm 



Apogee (Upper) Ring 


Naturally, the chute deployment is a condition that puts the rings £ ' 

while the landing phase puts both in compression. Roughly, both need to 
withstand joining the vehicle at lOg's in tension and 7.5 g s in compression, 
because me rings separate parts of the vehicle that have different masses, the 
forces each experiences will be different under the different conditions. 

The resulting applied stress on each ring will be given by 

aapp = F/A 


where F = ma 

where m = mass of the portion of the space- 
craft supported by the ring in ques 
tion under load 
a := acceleration in question 


The designer then compared this resulting applied stress to the critical stress for 
the case of a loaded truncated cone, which is given by 

2 . ,, i/2. 

a cr = (Y)(E)(t)(cosa)/(r)[3(l -v )] 


where 


Y = buckling safety factor — 0.33 
E = Young's modulus = 2.9 x 101 1 (Beryllium) 
t = thickness of cone = 0.005 m 
v = Poisson's ratio = 0.33 (non ferrous) 
r = smaller radius 
a = cone half angle 


The results for different loading conditions for the two Beryllium rings of Figure 
22 are summarized as follows: 


Ring 

r 

(m) 

a 

(deal 

Loading 

(a) 

Supported 
Mass (ka) 

Oapp 
(Mpo) _ 

Ocr 

(Mdo) 

Upper 

0.0475 

54.8 

7.5 comp, 
(landing) 

2.75 

0.28 

3550 

Lower 

0.100 

26.6 

7.5 comp, 
(landing) 

45.0 

2.17 

2620 

Upper 

0.0475 

54.8 

10.0 tens, 
(chute) 

87.0 

11.44 

3550 

Lower 

0.100 

26.6 

1 0.0 tens, 
(chute) 

34.0 

2.12 

2620 


Very obviously, the applied stresses in the rings are several orders of magnitude 
below the critical stresses, and therefore acceptable. 

Mass of rings 

Finally, we can roughly compute the masses of the rings, assuming that 
they are less in mass than a cylinder of similar radius and thickness. 

M = (density)(volume) < (p)(ft)(D)(t)(h) 

where for each ring 

p = 1826 kg/m3 
D = 35.0 cm 
t = 5.0 mm 
hupper = 9.0 cm 
hlower = 1 5.0 cm 


4 Summary 

This MSP structural design, then, represents a involved attempt to satisfy mission 
criteria for a spacecraft designed to return of a sample of Martian soil for 
Earthbound analysis. 

Although at points overdesigned, and in some areas not detailed enough, the 
structure (or at least its critical components) have been shown to be able to 
withstand some of the major loads placed upon them by the harsh and varied 
environments the craft would encounter on its round trip to and from the Red 
Planet. 
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Discussion of layout and propulsion basics, 


mission scenario 

Research on Viking, computer sketches 
Preliminary MSP sketches 
Assessment of MSP I, sketches for MSP II 
Computer drawings of MSP I, nozzle detail 
Finalized layout; components, shield, chute, 

gear . 

Discussion with Tim on component layout 
Layout of fuel, structure, fuel systems on 
computer 

Research on landing gear 
Sketches of gear, calculations 

Finalizing layout sketches 

rtkcArkor rnlrulntions 


More shock calc's 

Landing gear dan views, load calc s 

Final general layout details 

Decisions made on tasks finalizations or 
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Stress calc's on bus, gear, platform 
MSP final drawings started 
First text for final report written 
Final CAD work on MSP components 
Text writing for final report 
FirstLaser printing of report 
Pasted report 

Finished text, drawings of report 
Laser printed report 
Pasted, photocopied report 


TOTAL 61.50+ Hours for quarter 


'Involved group meetings (some not listed) 


Section B: 

Thermal Control and Guidance 


Terrance Yee 


SPACECRAFT THERM AI CONTROL 

Thermal design for this craft is primarily driven by the need to 
keep the propellent tanks between +10 and +55 degrees celcius and 
the batteries between 0 and +25 degrees celcius.(ref D6 ) The most 
critical conditions are during nightime on the surface and during 
ascent in the atmosphere. To remain within the allowable 
temperature range a variety of suface coatings, insulation, and 
active cooling and heating units are employed. 


Orbital Thermal Control 

Thermal control before descent is accomplished though proper 
surface coatings and orientation of the spacecraft, "ter release 
from the MOV and a short burn to precess in its orbit the MSP 
assumes a slow Votissary- mode of rotation about its axis of 
symmetry. The bottom of the hexagonal bus section is gold coate , 
the top of the hexagonal bus is painted white and the sides are 
vertically striped with 50% gold and 50% white paint. This paint 
scheme together with the shadow of the heat shield allow the 
spacecraft to alter the effective absorbtivity vs. emmissivity ratio 
by changing the spacecraft angle with respect to the sun. T e ang e 
of the symmetry axis with respect to the incident radiation is 
chosen by the onboard computer to give a daylight (non-eclipse) 
equilibrium temperature of 20C. If the spacecraft begins to warm 
the computer will tilt the top of the craft towards the sun there y 
shadowing some of the craft with the heat shield and exposing the 
white top of the hexagonal bus. This manuever will also be used 
prior to reentry to precool the craft to 15C which will lessen t e 
impact of the aerodynamic heating on reentry. If during orbit the 
spacecraft is too cool the top can be tilted away from the sun 
exposing the gold bottom of the hexagonal bus which will tend to 

warm the craft. 

This flexible effective surface was chosen due to the 
eccentricity of the martian solar orbit which alters the intensity of 
solar radiation depending on what time of year the mission is 



2 


performed during. It also allows us to be less precise in our 
calculation of the required mix of coatings. Those calculations were 


performed by analyzing the steady state radiation equilibrium 
setting the incident solar radiation and albedo radiation equal to the 
thermal radiation leaving the craft. Approximating the craft by a 
sphere of radius R the following results were obtained for the 


different mars-sun distances: 
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For an eclipse the following relations were used to find the 
time of the eclipse: 
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Using a lumped capacitance approximation for spheres of 
different radii the change in internal thermal energy is set equal to 
the energy radiated away during the eclipse. Specifying the 
allowable change in temperature as ten degrees celcius allows us to 
solve for the surface temperature. Once this is known a minimum 
ratio of insulation thickness to thermal conductivity can be derived 
from fourier's law of conduction. The actual thickness of insulation 
is seen to be greater than the minimum required so the craft will 
stay within the allowable temperature range. An analysis was also 
performed assuming no insulation and use of electric resistance 
heaters during eclipse as a comparison. 
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Reentry Heating 

The allowable temperature change on reentry is specified as 5 
degrees Celcius which allows a factor of saftey of 2 before the 
batteries overheat. Originally an ablative heat shield was 
envisioned but after considerable research and a more accurate 
determination of the reentry heating it was found that a radiative 
heat shield using a high temperature multilayer insulation would 
function quite well in this environment. 

The reentry heating of the shield was approximated by 
calculating the stagnation enthalpy and pressure at the stagnation 
point for numerous points during the descent then the heat transfer 
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rate was found using the chart provided. Setting the total heat 
transfer rate to the shield equal to the heat radiated away from the 
front surface allowed us to solve for the maximum temperature at 
this point on the front of the shield. This allowed us to choose the 
type and thickness of insulation to give the desired heat transfer 
and backface temperature. The beryllium honeycomb on the 
backface was not allowed to exceed 700 degrees Fahrenheit. It was 
also discovered that the thermal lag for a heat shield of about 45 
seconds allowed for an additional factor of saftey in the backface 
support structure since it will be at a lower temperature than 
equilibrium when it experiences the maximum drag force. (ref ) 

Two computer programs were written for this section. The 
first helped coallate the atmospheric, trajectory, and heating rate 
tables and charts. The second calculated a few stagnation 
temperatures outside the boundary layer in front of the shield 
assuming no dissociation (not a very good assumption) in order to 
get a rough idea of how much error is induced by assuming a 300 
degree Kelvin wall in using the given handout since a very high 
freestream temperature would make the differnce in heat shield 
temperature less important. (see appendix D1) 
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In the recirculation zone behind the heat shield the convective 
heat transfer is approximately one tenth of that experienced at the 
stagnation point. (ref 03 ) Approximating the temperature of the 
recirculation zone as 1400K and modeling the craft as a cylinder we 
get the following convective heat transfer relations for the bus 
during descent: 
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In order to maintain stability during descent a reinforced 
carbon- carbon "drogue chute" is deployed when entering the 
atmosphere. This cone- shaped body is tethered to the rear of the 
craft on either side of the motor by high temperature titanium 
cables. The drag cone is stored as the lid to the main chute canister 
and when deployed by coil spring will trail 10 meters behind the 
craft. Since the heat shield stagnation point was at 991 K we 
assumed this temperature as the worst case for the entire cable. To 
size the drag cone we assumed a 50 m/s gust which is 
approximately a 2% variation in drag force (about 100N). Letting 
this disturbance force act at 75cm from the center of the heatshield 
leads to a 75Nm disturbance torque. If we allow a maximum moment 
arm of the drag cone force of 1m then we require a minimum force 
of 75N which translates into an area of 220cm A 2 if we assume that 
the drag force per area of the heatshield is equal to the drag force 
per area of the drag cone. With a right cone (45 degree half angle) 
this translates into a surface area of 707cm A 2 and a base radius of 
15 cm. Adding a factor of safety to this, we increase the base radius 
to 21cm and the surface area to 980cm A 2. Making the cone out of 
5mm reinforced carbon-carbon (RCC) leads to a cone weight of 
784g. A simple stress analysis of the cable shows that a diameter 


of 3.175mm is capable of supporting more than 1000N force and 
weighs 360g. Therefore a total weight estimate of the stability 
system is 1.3kg allowing 156g for fixtures such as bolts and a 
spring. 


Thermal Control on the Surface 

After reentry the MSP may be warmer than desired, therefore a 
method of rapid cooling is required. The method chosen to 
accomplish this is three small louvers which are placed 120 degrees 
apart on the ascent body above the ascent fuel tanks. The louvers 
are opened by thermally set springs and held shut by latches which 
the computer opens upon landing and locks at night. They are top 
hinged, opening outward to expose an 8cm high, 3cm wide opening 
for the martian atmosphere to flow through. 
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Some concern was originally felt for the cooling effect of the 
expanding helium in the propulsion system. The expansion is solved 
using the ideal gas law: T.-z^k 
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This temperature decrease is certainly acceptable during the 
rocket firings while heat from the rocket motor and aerodynamic 
heating are offsetting the helium expansion effects. 

This brings us to the next problem area, the rocket motor. 

While the thrust chamber is ablatively cooled and has negligible heat 
transfer to the rest of the ship, the nozzle and throat are radiatively 
cooled and must be isolated from the rest of the ship. This is done 
using a three layered system consisting of a highly reflective outer 
coating of silver with a protective coating of SiO (150 Angstroms), 
a layer of beryllium as a heat sink, and a final layer of refrasil, a 
multilayer high temperature insulation. The analysis of this system 


follows: 
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The last ground problem is cooling at night due to radiation 
loss and convection. It is this area which requires the most 
insulation and indeed sizes the insulation used. A small UDMH 
decomposition fuel cell is used to warm the ship at night. 


S\«ir M*/S 


il\r $ ct~i e a + d n / ^ ej 1} t 


M - “f ,S x Jo 


H 




Cp - llo.lJ± 

■** > L K 


c \ ~ 


V 


f r r 7 T : ~ 0 , 552 


t re 5 ki + »on 


R< 


N, -- !< Tfr 


V^J 


^ / 

for Rp-nS 6 i 


^ z 
*'3 At 

OS’dl' } 

p I±r c 

Jj -y r„,, M 

( **i 

k *. ■■ 

,oos?y 


k 

r r/ 



„ ZO 1 ^ - 0.?^ 

r /V 3 6 / 


6 7 i ^ ‘* ,x 

> 1 


- . !«rj- 

, " " . 6 /^ 

frW 90^ 


or 

L 


k ^ J 
2 C \ 3 - l*° 


11 t US’ 

*% ’ 


T-UbV 


C£:r>, 


Wr 




ft 

TV 


3 ( 7 S 6 W 


2 ~ . i V U l.t^r^ 


M 2 .‘) y to ~ 3 ^- 


o v er 


9 


o \S°“ HI, til G? ~ l,3T7xlo ( 3 lt:ii 

g,Vf^ tUt ^ °' f <-o->W'« M li H(>xlo 6 j£^ -f or keroi^^ 

t'? esf-. «» »f> kyjrq^rwf' h-?s tvo«- s * i- ,, / + J,,\s valise 

J 

$ ‘ nr ^ t ^ ^ c too | vnlu<? tf v ft a v 4* Iq hi P* 

7~^ lends f-° 56 gr^^$ < 7 ^ Hy^?[n<r re^y*** 1 fc 

l^^p^cl^rr levels. 


c . ; : 


<? 


Ascent Heating 

Using the same methods as the descent the ascent trajectory 
was analyzed but with a nose made of RCC and .96cm of refrasil 
protecting the skin. Again as in the descent, the craft was allowed 
to cool to 15C before launch in order to alleviate heating 
constraints. 



Rendezvous Vehicle 


The rendevous vehical has very little on board that is 
temperature sensitive; the electronics are very simple and can be 
built to withstand wide temperature variations, there are no 
batteries or fuel, only superstructure, dirt, solar cells and an 
omnidirectional radio beacon. 

The surfaces not covered by solar cells are covered by gold to keep 
the temperature near 290K. 


fil JDANQE AND POWER SUPPLY 

Guidance is provided by inertial rate gyros and by sun and star 
sensors when outside the atmosphere. Final landing guidance is 
provided by a system which uses a prescanned map of the landing 
area which is stored in on board ROM computer memory. (see appendix 
D2) The map records safe landing sites and local landmarks for the 
entire probable region of landing. Just after the parachute is 
jettisoned and before the retro manuever a snapshot image is taken 
of the terrain and matched by landmarks and distinct features to the 
map stored in memory. Once the computer determines where it is, it 
can use the map to direct the craft to the closest known safe landing 
zone. A laser range finder is used to determine distances to the 
ground to allow the computer to make the decisions when to engage 
the main chute, when to jettison the chute, and when to engage and 
disengage the retro motor. 



The sun and star sensors along with inertial rate gyros and an 
accurate timer allow the ascent computer to determine the correct 
timing of the firing of the apogee kick motor as well as its correct 
orientation. Also the ascent computer will set a timer on the apogee 
motor that will pyrotechnically pierce the combustion chamber to 
effectively kill the thrust when the correct orbit has been reached. 

The on board computer actually consists of three parts. The 
first part is the descent computer which is the ROM board and the 
guidance for descent. This part is located in the hexagonal bus 
section and weighs about 1kg. The second part is the main or ascent 
computer. It is located in the ascent stage and weighs just under 
2kg. Its functions include powered flight guidance, directing the 
sample collector, and orbital equations. The third part is a subpart 
of the main computer. Its function is to periodically check all 
systems, especially thermal, during the dormant periods while the 
MSP waits for a particular orbital conjunction. It has limited 
control over the attitude jets and thermal system and will activate 
the main computer if a serious problem is detected. It has 
negligible power usage and negligible mass. 

Power is supplied to all but the rendevous stage by zinc-silver 
oxide batteries with a power density of 80 watt-hr/lb.(ref D2) 

There are two groups of batteries, one in the ascent stage and one in 
the hexagonal bus. The batteries for the sample collection motor are 
subdivided into many smaller cells in order to give a high voltage 
while the other batteries supply a lower voltage usable by the 
computers. The main computer system is assumed to operate at 
30watts and the descent computer at 8 watts, the monitor 
subcomputer is estimated to have such a low power usage as to be 
negligible. 

action 
ascent 

sample collection 
final descent 
deorbit and reentry 
adjusting precession 


power(W) 

time(min) 

energy(watt-hr) 

30 

35 

17.6 

330 

1 0 

55 

38 

5 

3.17 

30 

60 

30 

30 

40 

20 



burn to precession 
unforseen/reserves 


30 


1 0 


5 


1 2 


27.63 


158.4 

total 

Total battery weight is 0.9kg, 0.8kg is stored in the hexagonal 

The beacon on the rendevous stage is an omnidirectional 
beacon broadcasting a narrow band 100W signal ot 1 millisecond 
duration five times every second. To provide this power 145cm*2 
cross sectional area of solar cells must be used. On the 10cm 
diameter rendevous stage this translates into a 1 4.5cm strip around 
the cylinder of 11% efficient silicon cells at furthest distance from 
the sun If 17% efficient GaAs cells are used the strip will be 9.4 
cm wide. The remainder of the suface should be covered in gold to 
keep the craft warm. Silicon solar cells are currently available and 
their reliability is high so they were chosen for this mission, giving 
a weight of 0.15kg to the power system. 
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MSP THEFMA1 /POWFR/GUIDANCE/COOR DINATOR TIME SHEET 


•WEEK HOURS 

1 3 

2 5 

3 3.5 

4 5 

5 8 

6 1 0 

7 6 

8 1 5 


9 20 

10 25 


total 


100.5 
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boundary-layer SOLUTIONS 
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Keen found from a sin 

Heat-transfer rates and boundary-layer pro^^ #imilarity assump 

solution of the laminar hounaary-laye q £irat . orde r total differentia 

are used to reduce the equations to * ^solution is then found by an iterat 
equation, with variable coefficients. ™ * description of the solution is 

technique using 400-pomt profiles. A detaue 
also given elsewhere* 

, j * ♦v.i- '--maicarv -layer solution and > 

A number of assumptions are at^mlcVomposition is assumed to be 

be mentioned for completeness. - csscnce> thi8 means that multicom 
stant throughout the boundary ay * assumptions are exact for the 

diffusion is neglected. _ 1oC * 1 “™ rtic3 are continuously variable throu 

^:» -mi C al reactions are included by 

employing the reaction-conductivity concept. 

a « wall temocrature of 300 K* 

All calculations in the report have are rep orted here. Other cases 

Similarly, only stagnation-point so pressure gradients can be obt. 

more realistic wall temperatures and other pressu 

by application of this method. 
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Section C: 

Propulsion Systems 


Myles Baker 


Propulsion; 


1: Requirements: 

The Mars Surface Probe places several demands on its 
propulsion system. From the Mission Scenario, at least five 
separate impulses are required. 


Event: 


Initial orbital maneuvering 
Descent burn 
Landing retro burn 
Ascent from Mars surface 
Apoapsis injection 


Required Required 

Delta V: Thrust: 

<10 m/s < 10000 N 

<200 m/ s <10000N 

<200 m/s > 500 N 

—3920 m/s > 1000 N 

-v 830 m/s < 200 N 


In addition, the propulsion system should be as small and 
light as possible, since mass is at a premium in such a mission. 

It should be as simple as possible to increase reliability, and 
it should be very stable due to the approximately two years it 
will spend in flight to Mars. An additional restriction placed 
on the propulsion system is that it never cause the vehicle to 
accelerate at more than 100 m/s 2 , which is reflected in the 
maximum allowable thrusts. 

With these criteria in mind, the decision was to use two 
separate engines; one for the low-thrust apoapsis injection burn 
and the other for all other, higher thrust burns. Due to the 
varying requirements of the first four impulses, it was decided 
to use a throttleable, re-startable liquid propellant engine. On 
the other hand, because of the decision to have no active 
guidance on the final payload capsule, the advantages of a solid 
motor far outweigh its slightly lower specific impulse. These 


advantages include the lack of propellant control devices 
(valves, etc.), propellant tanks, and the fact that the pre- 
determined impulse delivered by a solid motor does not require 
any outside control mechanisms (other than a simple ignitor) . 


Liquid prop ellant Engine • _ 

The prime consideration in the design of this engine was its 
simplicity. There are several choices in the design of a rocket 
engine, each of which has a great impact on the final design of 
the vehicle. Some of the options, and selections for this 
mission, are listed below. 

Thrust: The primary restrictions placed on the thrust of 

the engine are the imposed design restriction of less than 100 
m/s* maximum acceleration and the need to overcome Mars' gravity 
with minimum gravity losses, which requires as high a thrust as 
possible. From Newton's second law of motion, 

p = m . a. 

* t ,max mm max 

With a burn-out ascent mass of 15 kg, a thrust of approximately 
1500 Newtons is required. 

Cooling: There are several methods of cooling a rocket 

engine, ranging from simple radiation cooling to regenerative 
cooling by pumping cool fuel through a cooling jacket around the 
engine. since this engine is rather small, the extra hardware 
and complexity required for regenerative cooling eliminate it as 
a viable option, while radiation cooling is very well suited for 
such small engines. However, due to the fact that the actual 
combustion chamber is located inside the vehicle during 


the mission, some additional cooling is required. Since there is 
a relatively short total burn time for the engine (approximately 
120 seconds) , an ablative combustion chamber liner will be used. 
Such materials are most reliable at chamber pressures of 1 MPa or 
lower, so 1 MPa is chosen as the maximum chamber pressure for the 

MSP liquid engine. 

Fuels: There are many available fuels for liquid rockets, 

and various parameters for selecting a fuel for a particular 
mission. in this case, our requirements are a high specific 
impulse, high density, and easy storeability for up to two years 
without maintenance. The latter requirement eliminates several 
of the higher-energy propellants, such as the cryogenics (Liquid 
Hydrogen) and the fluoridated oxidizers. On this mission it was 
decided to use Nitrogen Tetroxide as an oxidizer and 50% 
Unsymmetrical Dimethyl Hydrazine, 50% Hydrazine (50-UDMH) as a 
fuel. This propellant combination has several advantages, 
including hypergolic ignition (no separate ignition system is 
required), extremely high storeability, a relatively high 
specific impulse, and a broad base of engineering experience with 


such propellant combinations. 

Propellant feed: There are two types of propellant feed 

systems in use today: pump-fed and pressure-fed. The pump-fed 

system requires some additional machinery (a turbo pump, gas 
generator, and gas turbine) , but it allows lower propellant tank 
pressures, and hence thinner tank walls, lowering the total mass 
of a large vehicle. Pressure-fed systems, on the other hand, 
maintain the entire propellant supply at a high pressure 


and thus 


require heavier tanks, but virtually no machinery other than 
valves. For a vehicle as small as the Mars Surface Probe, the 
slight increase in tank mass is greatly compensated by the lack 
of supporting machinery, so a pressure-fed system is more 

appropriate . 

Materials: Recent advances in high-temperature, relatively 

oxidation-resistant composite materials makes the choice of 
materials for the motor structure fairly straightforward. Multi- 
cycle burn times well over the required 100 seconds have been 
demonstrated with engines constructed almost entirely from a 3- 
Directional Novoltex CR) type carbon fiber reinforcement m a 
silicon-carbide matrix 10 ' 11 . In the previously built engines, the 
fuels, flame temperatures, and motor sizes were similar to those 

in the Mars Surface Probe. 


solid Motor: 


The solid rocket motor for re-insertion into the Mars 
orbiting Vehicle's orbit is much simpler. It is only required to 
ignite, deliver a pre-determined impulse, and extinguish. The 
payload capsule will already be positioned in the correct 
attitude and spin-stabilized by the ascent stage, so no attitude 
control will be necessary. With spin-stabilization, if the spin 
rate is high enough, it is not necessary to adjust the thrust 
vector assuming a small uncertainty, so the engine mechanics can 
be extremely simple. The spin-stabilization of the uper stage 
will be provided by unwinding a pre-stressed ''clockspring" 
mechanism, in effect spinning the ascent and final stages against 


one another. 


The fuel used in the solid motor is also well-tested, off 

the shelf technology. It consists of Aluminum and Ammonium 

Perchlorate powders in a polyurethane binder. This compound was 

chosen primarily for its slow burning rate of 6 mm/s, which makes 

it possible to remain within the 100 m/s 2 design acceleration 

limit. To minimize the motor case mass and surface area, a 

spherical case will be used, with a nozzle very similar in its 

proportions to that used in the liquid engine. The materials 

used are also very similar to those used in the liquid engine, 

( R ) 

with the addition of a 4-Directional Carbon-Carbon Sepcarb 
throat insert to minimize erosion in this critical area. 


Propellant Tanks s_ 

There are two necessary sets of propellant tanks. The first 

set is used for the orbital maneuvering, de-orbit, and retro 
burns and is housed in the main bus of the vehicle, which is left 
behind on the planet surface. It contains approximately 5 liters 
each of N 2 0 4 and 50-UDMH pressurized to 1.2 MPa, allowing for 5< 
ullage, and 2 liters of helium gas at 25 MPa to pressurize the 
tanks. These tanks will be split up into three identical sets of 
three tanks each, with similar tanks (e.g. N 2 0 4 , 50-UDMH, and He) 
arranged symmetrically about the axis of the vehicle, 
small size of the tanks, and the relatively low pressure, 2014-T6 
Aluminum alloy (chosen for its good machining properties) will be 
used for the propellant tanks, and Ti 6A1-4V Titanium alloy (for 
its exceptionally high strength) for the Helium. Unfortunately, 


Nitrogen Tetroxide is not compatible with aluminum tanks, so an 
interior protective coating will have to be used. 

The second set of tanks, used only for the ascent burn, must 
contain about 19 liters of each of the propellants and 7 liters 
of helium, again allowing 5% ullage. The configuration for the 
ascent tanks is shown in figure Cl. with the materials similar to 
those in the descent tanks. The nested configuration was chosen 
in order to keep the ascent envelope as small as possible, and 
only weighs about fifty grams more than the optimum spherical 

tanks . 

Propulsion System An alyses: 

The design of a spacecraft is necessarily an iterative 
process, with the results from each part of the analysis 
affecting the parameters the remaining phases are based on. This 
necessitates a certain number of assumptions in the preliminary 
vehicle design, with succsessive refinements as the design 
progresses. The assumptions used in the analyses below have been 
through three such iterations, and so are fairly accurate. 


Liquid Propellant Engine; 

There are many parameters of interest in the design of a 
rocket engine, but perhaps the most important are the thrust 
(F t ), specific impulse (I sp ) , mass (m) , and dimensions, which will 

be developed here. 

The combustion characteristics of Nitrogen Tetroxide and 50- 
UDMH are tabulated in several references, each with slightly 



different values. 


The values below are 


taken from Sutton and 


Ross 


CC - 21 


3100 K 
396 J/(kg-K) 
1.24 

1450 kg/m 
910 kg/m 

1.62 

1.01 


Flame Temperature: 

Gas constant (R) of Exhaust: 

Ratio of Specific Heats (Y) : 

Density: N 2 0 4 : 

50-UDMH: 

Mixture Ratio (0:F) 
by Mass: 
by Volume: 

Table X.X: . . . _ . 

Combustion Properties of Liquid Fuels 

Nozzle: Assuming a ratio of €=20 between the nozzle exit 

area (A e ) and throat area (A t ) , the mach number (obtained through 
a computer iteration) at the exit plane is found to be M e =3.93. 
Given the chamber pressure (P c =lMPa) , the pressure at the nozzle 
exit plane (P e ) can be obtained from the equations for isentropic 

flow: 

P. - P e X [1 ♦ 

p e = 10 6 X [1 +li^-(3.93) ! f W< '"'’ > 
p e = 4439 Pa 

Given the chamber pressure, the exit pressure, and the back 
pressure (P b ) , the thrust coefficient, defined as 

F t = P c A t C f 


-v- ^r " 1 


can be obtained from the relation 

which , for this engine, at the mean surface pressure at the MSP 

landing site of 700 Pa, is 
C 


"f , ideal 


f ZluM 
|. u^\ 


J a m 

* . _ . \ 1 1.14-n 

z -\ ■ m-i 

\ 1H H j 

I ’ \itCCCCC ) 


i CoC GlC 


(■ 2 . 6 ) 


C, - A , = 1.77 

, ideal 


Assuming a bell nozzle profile of maximum efficiency (a perfect 
bell nozzle of expansion ratio 30:1 truncated to 20:1), the 
nozzle losses (due to friction and nozzle divergence) are 
approximately 2%“ \ so the actual thrust coefficient is 
C f - X 0.98 - 1.73 

Given a maximum thrust of 1500 N and a maximum P c of 1 MPa, the 

required nozzle throat area is 

A t = F t /(P c C f ) = 1500/(1.73 X 10 6 ) 

A t = 9 cm 2 

For a throat area of 9 cm 2 , the inside diameter at the throat is 
3.4 cm, the exit diameter is 15.1 cm, and from the bell nozzle 
optimization curves [C ' 31 , the nozzle length is 18.6 cm. 

The mass flow through the nozzle can be obtained by a mass 
balance at the throat, where necessarily sonic, so 

• x m / Y I 2- Vor-0 

m “ A t P cV^T t l,YM ) i i±±L 

a -4 ?. / i it- 2 ( 2 Vo-»-o 

m = (10 Pa) (10 m )J(vcb\i X (fi)\ VLi M ) 

m = 0.533 kg/s 


at a maximum thrust of 

F = C f A t P = (1-73) (10‘ 4 ) (10 6 ) 

1 t ,max f t c 

F = 1560 N (at surface) 
t,max 1575 N (in vacuum) 

This gives the liquid engine a specific impulse of 

■*-sp a, 

Isp = 298 sec (at surface) 

301 sec (in vacuum) 

Combustion Chamber Geometry: In order to minimize pressure 

losses in the combustion chamber, the velocity just upstream of 
the throat should be very small, implying a large chamber to 


throat area ratio (A c /A t ) . In this case, A c /A t will be set to 
6.0. Thus the radius of the combustion chamber is 
r c = /6A c /;r 

r c = Vo. 0054/*' = 4.45 cm 

A second parameter in the design of a liquid-fuel combustion 
chamber is the characteristic length (1*) , which is defined as 
the chamber volume divided by the throat area. This length is a 
useful measure of the amount of time a given element of the fuel- 
oxidizer mixture remains in the chamber, and must be high enough 
to ensure complete combustion. In the reaction control thrusters 
used on many early American launch vehicles, Monomethyl Hydrazine 
and N A were used as propellants and the characteristic length 

[C-2] This thrust chamber was used as a model, 

was approximately 0 . 5 meters 

requiring a combustion chamber length (approximating the chamber 
as a cylinder) for the MSP engine of 

l c - l*(A t /A c ) 

1 = 0.5/6 = 8.33 cm 

Allowing for 5 cm clearance at the injector for fuel lines and 
the thrust vector control (gimballing) mechanism, the total 
liquid engine dimensions are: 

D max = 16 Cm 

Itot = 30 CItl 

Ablative Cooling Insert: Ablative materials are often used 

in cooling small, short duration thrust chambers. When exposed 
to the heat of combustion, the ablator, a fiber reinforcement in 
an organic matrix, decomposes, leaving a porous, insulating layer 

The decomposition of the organic matrix 


of char on the surface. 


is an endothermic process, and the pyrolysis gases released tend 
to insulate the remaining material from the hottest of the 
combustion chamber gases. 



Figure C *4 : Ablative Chamber Cooling 

A silica fiber reinforcement in a phenolic matrix will be used as 
the ablative liner of the combustion chamber. The depth of the 
porous char layer, and hence the necessary thickness of the 

ablator, is given by the relation cc 31 

d c = /burn time in seconds mm 
d = J 120 1 mm = 11 mm 

C 

This analysis does not consider the fact that the engine will 
reach ambient (spacecraft) temperature between firings, so the 
actual char depth should be significantly less. Nevertheless, 
adding a margin of safety, the thickness of the ablative material 

should be approximately 1.5 cm. 

Propellant Injectors: Given the overall mass flow rate, the 

mixture mass ratio, and the propellant densities, it is easily 
seen that the necessary volumetric flows for maximum thrust are 
Q f = 0.225 1/s (50-UDMH) 

Q o = 0.227 1/s (N 2 0 4 ) 

and, given an orifice-type injector, the volumetric flow rate is 


given by the equation 


[C-4] 


1 1£> F7 


Q - A ,cfi ^ 

where C, is an experimentally determined discharge coefficient, 
usually between 0.6 and 0.9 (assume 0.75), and P, is the 
pressure drop across the injector. For a set of n injectors, the 

total area is 


U V2A?; 

which is, for a 200 kPa pressure drop, 


A- , = 1.43 X 10' 5 m 2 

i.f 

A- = 1.82 X 10‘ 5 m 2 
1 r 0 

For efficient combustion, it is desireable to have as many 
doublet-impinging pairs of injectors as possible in order to 
atomize the propellants effectively. This is limited, however, 
by the machining of the injector manifold and the injectors 
themselves. One possibility is a set of 19 injector pairs 
arranged as shown in figure Q.l. The diameter of an individual 

injector is 


D : 




A 7T 


so, for 19 injectors, the diameters for the fuel and oxidizer 

injectors, respectively, are 

D- „ = 0.979 mm 

D- = 1.104 mm 

1 ,0 

The injectors should be positioned at angles such that the net 
momentum of the fuel and oxidizer streams is parallel to the axis 

of the engine, so 

O-'-'-pr 

. -1 r 

= Sin _d_ 

LA:'' 


fo 


0 v o $in(<p o ) m,v f sin(^ f ) 

j rpujfc r>o A; c ,.o 

1 T , ^ j 

J 


1 




Assuming a fuel injector angle of 20° from the engine axis, the 
oxidizer injector angle should be 

c . -i r i o^x ( mi) ( 

= Sin /( ^lo ){ M'S / J 


?° 

y° 


= 15. 5 C 


Chamber Stress Analysis: The stress analysis of the 

combustion chamber is very straightforward, considering only the 
hoop stresses. For 3-Directional Carbon-Silicon Carbide 
Novoltex (R \ p= 2300 kg/m 3 , £T y = SOMPa, and for hoop stresses, 
t H = r c (P c /c y ) X S.F. 

where S.F. is the safety factor (1.5) and r c contains the 
ablative chamber liner. Substituting in the values, 

t = (0.0445 + 0. 015) ( 10 6 / (80 X 10 ) ) X 1.5 

W V 

t = 1,1 mm 

w 

A combustion chamber with walls this thin would be extremely 
difficult to construct reliably, so an additional safety factor 
will be added to the wall thickness by increasing it to a uniform 

2 . 0 mm . 

Mass Estimates: Given the material properties and 

dimensions calculated above, it is very straightforward to 
compute the mass of the engine as the sum of the masses of the 
components. The volume of 3-D Novoltex"’ in the combustion 
chamber and nozzle shells is 


V 


c/c 


= (0.002 m) ( 2 ?T ( r c + 0.015 m) + A n ) 


Where A n , the nozzle surface area, is given by 
A = 70 X A t 


[C - 21 


V = ( 0 . 002 ) (2^(0.0445+0.015) + 70(0.0009)) m 

C/C 


V 


c/c 


8.74 X 10‘ 4 m 3 



/ (Itaronce 


so the mass of the chamber and nozzle shells is approximately 
M . . . = 2300 kg/m 3 X 8.74 X10 4 m 3 

1A shell 

Mshell = 0.20 kg 

The ablative insert mass, assuming a uniform density of 
kg/m 3 is 

M abl = 2 *< r c + 0 • 015/2 ) X l c X abl X 0.015 

Mabl = 2 tl ( 0 . 0445 + 0.0075) X 0.0833 X 1740 X 0.015 

M abl = 0.71 kg 

And the mass of the injector plate, assumed stainless steel (p- 
7600 kg/m 3 ), 1 cm thick, and 50% hollow for fuel manifolds, 

M =7T(r + 0 . 015 ) 2 X 0.01 X 7600 X 50% 

ll inj ' c 

M. ■ = 0.80 kg 
A i n j 

Thus the entire engine mass is approximately 

M = 0.20 + 0.71 + 0.80 = 1.71 kg, 

eng, l iq 

Allowing some excess for throat reinforcement, plumbing, and the 
gimballing bearing, 


sol id orbital Insertion Moto r : 

The analysis of the solid motor is very similar to that for 

the liquid motor. The properties of the solid Polyurethane/ 
Ammonium Perchlorate/Aluminum fuel, again from Sutton and Ross [C 


2 ] 


are : 


Flame Temperature: 

R of Exhaust Products: 
Y of Exhaust Products: 
Burning Rate: 

Density: 


3088°K 

284 J/ (kg-K) 
1.17 

6 mm/s 
1750 kg/m 3 


Dimensions: Assuming a chamber pressure 


of 1 MPa and a 


nozzle exit area ratio of 20:1, the Mach number at the exit plane 
is found (again by computer iteration) to be 3.65. At the exit, 
P. _ ( 10 ‘)[1 + !dil3.65 ;: ) 

p = 5452 Pa 
e 

And the thrust coefficie nt i g 

, ideal 


Cienwa ~ 

hAuifl l \«r t -/g£LV- ,T 1 

~'\j i.n ” Vi-n H / l [loooooOf \ 


&ri- 0 (zo) 

(OoOOCO 


C- . , , = 1.849 

W, ideal 


Again applying the 98% nozzle efficiency (similar geometry to 
liquid engine) : 

C f = 1.812 

Assuming a maxrmum thrust of 200 N (2.0 Kg final mass at 100 


m/s 2 ) , 


200 tJ 


j.?U ( 10 * 


At “ C (^ 

A t = 1.103 X 10' 4 m 2 


A t = 1 cm2 


So the actual maximum thrust is 

F 


’ t ,max 


= (10' 4 ) (10 6 ) (1.812) = 181.2 N 


The mass flow is again given by 


m = a A p = 

r -a f ) f 2 

m = (10 ) (10 (i a <-1 i 


m = 0.0687 kg/s 

And hence the specific impulse. 


I = 181. 2/( (0.0687) (9.81) ) = 269 sec 


ideal rocket equation, gives the 


The specific impulse, from the 
propellant mass and volume. 


M p = 2.0 kg X [exp 


M p = 2 [exp 




‘i* 


- i] 


i] 


M = 0.734 kg 

p 

v P = v<?p = 4,19 x 10 4 

The burning area (A,) tor these conditions can be found from 

A b,max = */((\t b ) 

A = (0.0687)/((1750) (0.0006)) 

b , max 

- 65 -4 cm ! 

orain Design: For a neutral-burning, cylinder-based grain 

configuration (such as the conventional star grain), the cylinder 
length and radius, in order to meet the criteria on area and 

volume, must be 

r = 2 Vp/A b = 2(4.19 X 10 4 ) /6 . 54 X 10 3 

r = 17.1 cm 
1 = Vp/(/rr 2 ) = 6 mm 

obviously this grain configuration is absurd, so another must be 
sought. Since this motor operates in free fall, a neutral burn 
is not required, and the options for grain configurations is 
almost endless. This motor will have a "sphere within a sphere" 
grain configuration as shown in figure C-4- ™e burning area as 

a function of time is 

A jV t l - (r„ -O j 

and the remaining propellant mass is 




\ 

J 


Using the fact that the thrust in vacuum is directly proportional 


to Ah and P , a computer spreadsheet was programmed to calculate 
the chamber pressure and acceleration as functions of time 
(figures x.x and x.x), and various dimensions were tested. The 
final propellant grain dimensions are a total gram radius of 
4.65 cm and an initial bore radius of 1.0 cm. As can be seen 
from the graphs, the maximum acceleration is 90.0 m/s 2 and the 
maximum chamber pressure is 1.25 MPa. This pressure is slightly 
higher than initially assumed, but easily feasible. 

Stress Analysis and Mass Estimate: The stress analysis is 

nearly a duplicate of that used for the liquid engine except for 
the fact that this case is spherical rather than cylindrical. 

Chamber wall thickness is given by 

t H = r c /2 X (P c /<T y ) X S.F. 

With similar materials, 

t u = (0.0465/2) (1.25/80) (1-5) 

t = 0.545 mm 

W 

This chamber is also impossible to reliably construct, so g ' , 

a 2 mm uniform wall thickness is assumed. This gives the total 

mass of the case and nozzle . 

M = ( 0 . 002 ) (4 7i( 0 . 0465) 2 + 70 ( 0 . 000 1 ) ) ( 2 3 00 ) 
tot ' 

M tot = 0.157 kg 

Allowing 0.093 kg for throat inserts, internal insulation, and 
attachment hardware, the mass of the solid motor is 

^fueled = O' 984 k< 3 

M dry = °‘ 25 ° Kg 

Solid Motor Spin-up Mechanism (Clockspring) . In order 
have stability in a spin-stabilized spacecraft, the rotational 
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kinetic energy must be much greater than the maximum anticipated 
disturbing went. Assuming a 1 mm misalignment of the thrust 
vector with the center of mass at 200 H thrust, the disturbance 

torque is 

^ = (0.001 m) (200 N) = 0.2 N-m 
The rotational kinetic energy of a spinning body is given by 

E K = I <*> Z /2 

where I is the moment of inertia about the spin axis and ul is the 
spin rate in radians per second. Assuming the rendezvous stage 
to be a uniform cylinder of r - 5 cm, m - 2.7 kg, then I is given 

by 

I = mr 2 /2 = (2.7 kg) (0.05 m) z /2 
i = o. 00375 kg - ™ 2 

Assuming a kinetic energy of five times the disturbance torque, 

/iCKO.Z. M-rn) 

w — y o oo^r^-n*' 

03 = 23.09 radians/sec 

With a safety factor of over 2 to include the depletion of 
propellant mass, 

= 50 radians/sec 

The angular acceleration of a rotational mass is 

r . 

« T" ’ 

and the torque exerted by a rotational spring is 
T = -k & 


+ 1 e = o 

dtS l 

& = -a 0 cos(Jf t) 

^ = rz - sin Wr } 


which implies 


At the maximum 


rotational velocity, this is 

k)max = 

K I 

Assuming a spin-up angle of 10 revolutions (62.8 radians), 

* = 

k = 0.00238 N-m/rad 

This is a very small torque, and a spring to provide it, along 
with a bearing and shaft, should weigh no more than 0.1 kg. 


Liquid Propellan t Tanks:, 


Descent: During the descent phase of the mission, there are 

three separate liquid motor firings: orbital Maneuvering, De- 

orbit, and Retro Landing. The mass of fuel required for each 
burn is calculated by applying the ideal rocket equation, 


M 


fuel 


= iw e *p(-$vV 11 


using the v's and M >rt 's at the given time. Applying the 
equation , 


M 


f , om 


0.802 kg 


M' = 7.353 kg 
M/ ° = 1.850 kg (max) 

f , ret 

H (J =10. 105 kg 

Ll f ,desc 

Allowing almost 2 kg of fuel for ullage and contingencies such as 
an extended retro firing, the total descent fuel mass is 12 kg. 
Using the mixture ratio and densities of the propellants, the 
necessary tankage is 5.1 liters each of H A and 50-UDMH. The 
volume of Helium necessary to keep these tanks pressurized 

be calculated by modelling the process as an adiabatic 


MPa can 


expansion of an ideal gas where 


PV^ = Constant 

For an adiabatic expansion from P, , V, to P,. v 2 - the required V, 

V, - V, X p/VlP^- P,*> 
v = (2 X 5.1) x 1.2” 1 -“/<25''’-“-1.2 ,/, -“> 

HC 

V u =2.0 liters 

He 

Taking three identical sets of tanks, the radii are 

r f = 7.4 cm 

r = 7 . 4 cm 
0 

r He = 5.4 cm 

Using 2014-T6 Aluminum ( p = 2800 kg/m , cT y - 386 MPa) for 
propellant tanks and Ti-6A1-4V Titanium ( p= 4430 kg/m 3 , <T y = 
999MPa ) for the Helium, the required wall thicknesses are, from 
an analysis identical to that used in the spherical solid rocket 

case: 

t f = 0.17 mm 
t = 0 . 17 mm 

o 

t He = 1.01 mm 

The thicknesses of the fuel and oxidizer tanks are thinner than 

can be practically manufactured, so assuming a minimum 

. ^ o Trim this total tank mass is 

machineability of aluminum to be 0.3 mm, 

M t = 12 TT ( r f 2 t f ^> Al + r 0 t o p Al + ^He t HePT i ' 

M = 12 7T ( (2) (2800) (0.074 2 ) (0.0003) 

1 V (4430) (0.054)^(0.00101)) 


M t = 0.84 kg 


Adding over 50* for valves and plumbing, the total fuel system 

mass for the descent stage is M = 1.50 kg. 

Ascent: The fuel requirements for the ascent burn cannot be 

calculated from the ideal rocket equation due to the effects of 
aerodynamic drag, but since the burn time and thrust are known, 
the mass of fuel expended can be calculated using the specific 


impulse : 


M P cjo T-.P 


Where t b is the burn time in seconds. For the ascent stage, with 
a 1500 Newton constant thrust for 82 seconds, 

M = (1500) (82)/ ( (9.81) (298) ) 

P 

M = 42.0 kg 
p 


Adding 5% for ullage, 


M = 44.1 kg 

using a breakdown similar to the descent tanks, this requires 
volumes of 

V f = 18.5 liters 
V =18.9 liters 

O 


V 


He 


= 7.3 liters 


in order to minimize the envelope of the ascent vehicle, the 

tank configuration shown in figure C.7 was selected. The 

spherical caps and the cylindrical section have a radius of 

R =?/!£ 

K V W 

R = 18.5 cm 

and the cylindrical section (containing N,0 t ) has a length of 

L = vy&rR 2 ) 

L = 18.0 cm 





The Helium tank has radius 

_ Z (o.oori)' 

“ V 


He 


R He = 12.0 cm 


Again using 2014-T6 Aluminum for propellant tanks and T1-6A1-4V 
Titanium for the Helium, the required thicknesses are: 


t , = 0*43 mm 

^sphere 


"*cyl 


t H e = 


0.86 mm 
2.26 mm 

The total mass of the ascent tank structure is 

M = -7r((6R 2 t sphere + 2RLt cyl ) (p Al ) + 4R He 2 t He p Ti ) 

M - nr( 6(0 185) 2 (0. 00043)+2 (0. 185) (0.180) (0.00086) ) 
*£ (2800) *ntW (.120)^(0.00226) (4430) 

M = 3.07 kg 

, . ■ ,rai voq the final tank mass is 

Adding 0.43 kg for plumbing and valves, 


M = 3.50 kg 
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Section D: 

Orbital Mechanics 


Bob Langberg 


Orbital Trahsf 


AMP 


Aacmri/niiscBiiT Trajectories 


Charac t i^istics Pebtikfht to Orbxt^x . 

6.4182 10 23 (kg) 

universal Gravitational constant 6.6732 10" 11 W/*<3 s2) 
„ Gravitational Constant tor Mars ^“(iSs/saf “ 

Equatorial Gravitational Acceleration 3.73 (m/s2) 

Mean Radius 
Sidereal Day 


Angular Velocity 

Rotational Speed (Equatorial) 

(28.5° Latitude) 


3402 (km) 

24 hours 37 minutes 23 seconds 
24.623 (hours) 

1477.4 (minutes) 

88643. (seconds) 

0.24367 (deg/min) 

4.2529 10" 3 (rad/min) 
14.468 (m/min) 
0.24113 (m/s) 

12.715 (m/min) 
0.21191 (m/s) 


Q»BITXHG VeHX C T. E fMOVl OrBITM. P^lAHETl^ S 


Circular 

Ascending Node at 250 Longitude 
Inclination 30° 

Altitude 600 (km) above mean radius 
Orbital Radius 4002 (km) 

Orbital Velocity 3.2715 (km/s) 

Orbital Period 2 hours 8 minutes 6 seconds 

2.1351 (hours) 

128.11 (minutes) 

7686.4 (seconds) 


T.XMDIMG SxTE 


Depression Hellas, 
Altitude 
Radius of Site 
Longitude 
Latitude 


near Crater 29 

4.2 (km) Below Mean Radius 

3397.8 (km) 

283 (Deg) 

-28.5 (Deg) 
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initially the descent trajectory was modeled without 
taking the atmosphere into account, using ideal eguations 
tor elliptical orbital motion (Kepler Motion). In this 
analysis there was an orbital plane change, calculated with 
spherical trigonometry, at the same time as the deorbit 
burn. This combined burn put the Mars surface Probe (MSP) 
on an elliptical orbit that intersects the Martian Surface 
at the landing site. It was optimized with the constraint 
that the flight angle at landing could not be less (not 
shallower) than 15’. Orbits of eccentricity varying from 
0.08 (nearly circular) to 0.93 (an extremely narrow ellipse) 
were examined. The delta v retired was least at the more 
circular ellipses, with an absolute minimum occurring at the 
ellipse with eccentricity of 0 . 14 , which had the minimum 
plane change angle (10.1'), unfortunately its angle with the 
horizon upon landing was 7 . 83 ' which was too shallow. 

The best case which did meet the 15* requirement was 
the transfer ellipse with eccentricity of 0.270 . It had a 
velocity at apoapse (4002 km) of 2.795 km/sec, which is 
0.476 km/sec less than the MOV's circular orbital velocity 
of 3.271 km/sec. This orbit takes the MSP through a 58.7 
change in true anomaly, to land with an impact velocity 
3.409 km/sec. The entire deorbit takes 21.8 minutes, during 
which time Mars rotates 5.3l‘ on its axis. The plane change 
angle for this case was 11. 7’, which when combined with its 
velocity change for deorbit (0.476 km/sec) required a total 
delta v of 0.781 km/sec. This was the best that could be 
done with the plane change while neglecting the atmosphere. 


Avoiding the Plane Change Maneuver 

Jter the ideal analysis was completed' it was decided 

that a significant fuel savings (305 m/sec) could be 
attained by avoiding the 11. 7" plane change. By eaving the 
MOV orbit with a small delta v of 8.1 m/s (0.2% of e 
circular velocity) the MSP would advance or recede y a ra 
of 2 S’ per orbit. After 72 orbits, the MSP may move up to 
+ 180 ’ , allowing it to be positioned anywhere on the orbit 
With precise positioning it may then deorbit as Mars rotates 
and the target point passes below, without the necessi y o 
a plane change. This positioning would take up to 6.243 
Mars sidereal days <6 sidereal days and 9 hours) and would 
allow the MSP to position itself at any point in orbit to 
successfully deorbit and land at the specified sampling 
site. The extended time in orbit presents no apparent 
thermal problems. 

Limits of Atmospheric Effects 

To accurately model the effects of the atmosphere, th 
deorbit trajectory was assumed to be a Kepler ellipse un i 
an altitude of 250 km was attained. The drag force on the 
MSP at an altitude of 250 km is 1.04 lO" 4 N. The maximum 
drag is 6.92 10+ 3 N which occurs at an altitude of 30.1 
This difference, in excess of seven orders of magnitude 
shows that not only is it a very reasonable assumption to 
ignore atmospheric effects above 250 km, but that 
atmospheric drag probably could have been ignored until 
about 150 km altitude where the drag is approximately 0.6 N. 

Incremental Time Step Analysis 

kt 250 km, the velocity and flight path angle was 
computed for the ideal transfer orbit. These were used as 
initial conditions for a time stepping analysis. In thi 
analysis, the aerodynamic force on the MSP, which was 
assumed to be drag only, with no lift is computed. The drag 


force was computed based on a Cd of 1.0 which is referenced 
to the projected area of the MSP. The gravitational 
acceleration was computed as a position of radius us '-" g 
standard inverse square law at each position during e 
descent. These two forces were used to compute the 
acceleration in the radial and tangential directions and the 
resultant changes in velocity and position. Time mere 
ranged from one second to ten seconds, with five secon 
increments used during the period of maximum aerodynamic 
forces. The MSP mass was changed as the heat shield is 
dropped as the parachute is deployed. The Cd reference area 
was changed at the time of parachute deployment, but the cd 
is kept at 1.0. Both of these changes take place 
approximately 1.2 km above the landing site. 

Terminal velocity of Heat Shield and Parachute 

Terminal velocity design charts were generated 
heat shield and parachute design analysis. These charts 
were used in the sizing process only and were not directly 
part of the incremental time/dynamic analysis. The fina 
dimensions for the heat shield and parachute diameters are 

1.5 and 4.0 m respectively. 

Atmospheric Model . 

The atmospheric model used to calculate the vanati 

of density with altitude, used for drag calculations is 
based on: The Mars Reference Atmosphere , 1982, pu l 

The committee on Space Research, Chapter 1: "Post-Viking 

Models for the structure of the summer Atmosphere Mar , 
written by A. seiff. Data from this source concerning 
atmosphere in the southern hemisphere, during summ , 
altitudes was used to construct an exponentially d » oa ^ 
atmosphere model. The model is based on the density a 
mean surface, 1.78 10’* kg/m3 , and an exponential l 
height, 11.75 km, where the density has fallen off by 

density(alt) = density(O) * expC-alt/scale ht.] 


This model is most accurate where it has the greatest 
influence on the HSP, at low altitudes where the density is 
greatest, and the resulting aerodynamic forces are largest, 
see graphs for the shape of the atmospheric models. It is 
significant to note that the density has fallen off e 

linear scale at 70 Km altitude. (0.26% of the density at 

surface. 

optimizing Deorbit Delta V with respect to Eccentricity 

The transfer ellipse that sets the initial conditions 
for the descent analysis was originally assumed to be the 
transfer ellipse from the ideal case with an eccentricity 
0.27 . Due to the drag on the heat shield and parachute, 
the landing trajectory is nearly vertical near the surface. 
Thus the IS" restriction on landing angle is no longer a 

constraint for the real descent case. 

By doing a computation comparing the deorbit fue 
required vs. the eccentricity of the exoatmospheric transfer 
ellipse , a simple relationship was developed. It shows that 
the delta v only depends on the original circular velocity 
and the transfer ellipse eccentricity, e. 

Delta V = V Circular * [ 1 - Sqrt(l-e) ] 

This function is tabulated and graphed. It was found that 
another significant fuel savings could be realized ^mating 
an even shallower approach to the Martian surface, and t 

shallower the approach, the more fuel saved. 

The problem with making the approach more shallow is 
that the distance traversed becomes greater, and the 
accuracy of the reentry is reduced. For this reason the 
transfer ellipse was made only slightly more circu ar. 
original eccentricity was 0.27, requiring a delta v of 0 4 
Km/sec, and the new eccentricity is 0.10, requiring a delta 
v of only 0.169 km/sec. This causes the true anomaly, 
traversed in descent to increase from 53.42 to 96.10 , an 



the time of descent to increase from 22.37 to 36.04 minutes. 
However, most of this increased travel is exoatmospheric. 

The atmospheric true anomaly traversed (250 km to surface) 
increases from 11. 23* to 16.00”, and the time in the 
atmosphere (250 km to surface) goes from 5.85 to 7.19 

minutes. 

These two trajectories are essentially the same near 
the surface where the velocities and flight path angles are: 
161.7 vs. 159.9 m/sec and 66.8” vs. 68.8”, for the old an 
new cases respectively. These values are for an altitude 5 
km above the surface. Thus, for this analysis they are 
essentially compatible, and since the e = 0.10 case is more 
fuel efficient, (0.169 vs. 0.476 km/sec, a 64% savings) i 
was chosen as the final descent trajectory. The fuel 
savings was approximately 10 kg. 

Range Error Due to Uncertainty in Deorbit Delta V 

Because we enter the atmosphere at a rather shallow 
angle, (17.6” at 100 km altitude), the accuracy of the 
trajectory was questionable. Since the vehicle is in free 
flight with only small attitude control thrusters for 
orbital corrections, it is of importance to know how 
sensitive the trajectory is to the only parameter we can 
effectively control, the deorbit delta v. The same time 
stepping analysis was run for a 2» decrease in delta v. The 
new case placed the MSP 13.64 km short of its ^ginal 
target, but this is a small change considering the fact th 
the deorbit trajectory traverses over 5800 km. 



The final descent trajectory has the following 


characteristics : 

Eccentricity 

True Anomaly Traversed 


0.10 

( exoatmospheric/atmospher 
82.1° /16 . 0 


. . 4002 km 

Radius at Apoap 3.102 km/sec 

V Terminal <300 m above surface) 

Delta V (combined deorbit & retrofire) 


Max. 

Max. 

Max. 

Time 


«ach (100 km alt ) 19-3 

Dyn. Press. (30.1 km air.; ' 

Deceleration (30.1 km) 8*2 earth g 
during max deceleration (3000/250 m/s) 

/ D S GC 


^WCEHT ^HAT.YaXfl 


The ascent was modeled the same way as the descent. 

The assumed Cd was again 1.0 . The gravitational 
acceleration and the atmospheric density varied with radial 
position. Because the 82 second 1500 N ascent burn cannot 
be assumed to be impulsive for a ten minute flight, the 
thrust was considered to act in a constant manner and was 
just considered to be another force on the body for the 
duration of the burn. It was assumed that the mass 
decreased linearly during this period. The resulting 
acceleration was computed each second during the burn and 
every five or ten seconds after burn out. These 
accelerations were used to calculate the new velocity and 
position in polar coordinates, just as they were in the 
descent analysis. The basic parameters that were varied 
were the length of the 1500 N burn and the launch angle. 
From these initial conditions, the delta v required at 600 
km to circularize the orbit at the same velocity as the MOV 

was computed . 


Optimum Ascent Angle and Minimum Insertion Delta V 

The required delta v was quite sensitive to laun 
angle, as measured from the horizon. The optimum angle 
varied between 60° and 62° in a nonlinear fashion, depending 
on the duration of the ascent burn. At each burn duratio 
the optimum angle was determined to within 0.5 . By 
launching at 0.5° too low an angle, the ascent stage fell 
approximately 25 km short of its 600 km target altitude, 
launched at 0.5° too great an angle, it would reach 600 km 
with a flight path angle of about 10°, which must be 
corrected for with a larger delta v for circularization. 

The optimum launch angle allows the ascent stage to Dust 
reach the 600 km mark with a flight path angle as close to 
zero as possible to minimize the delta v required. 



It was also found that with longer launch burns the 
velocity at 600 km was greater (at the optimum launch 
angle) , requiring less and less insertion/circulanzation 
delta v. The penalty here is the increased launch mass for 
fuel and tanks. There is obviously an optimum distribution 
of the increased mass between the boost stage and the 
apoapse stage, but we had a driving design condition other 

that total mass to consider. 

The driving criteria for the launch trajectory urn 

out to be the apoapse kick solid rocket motor. With its 
burn characteristics determined by the spherical gram 
configuration, the maximum delta v it could provi e wi 
exceeding 10 g's was between 800 and 850 m/s. With our 
original launch mass of 46.4 kg the minimum possible de 
at the optimum launch angle was 1481. m/s, which could not 
safely be considered. A number of alternate cases were 

investigated. 

The final optimum launch profile increased the launch 
mass to 56.8 kg (a 10.4 kg increase), neglecting the 
increase in fuel tank mass. This profile required a launch 
angle of 61.0- and a burn for 82 sec of 1500 N. The ascent 
stage reaches 600 km altitude after 10.8 -notes with a 
flight path angle of only 2.92'. The required delta 

828.0 m/sec. . , 

If the flight path angle had been 61.5 the required 

delta v would have been 930.4 m/sec (a 102.4 m/s increase) 

and if launched at 60.5’ the ascent stage would have fallen 

about 40 km short of its 600 km target. 


Rendezvous Strategy 

For rendezvous purposes the ascent stage cou 
placed into an orbit with a slightly different period than 
that of the MOV, thus letting the sample return ve lc e 
catch up to the MOV without the MOV having the burn any 
additional fuel (passive recovery). As discusse or 
deorbit positioning, a ±8.1 m/sec delta v would let 


ascent stage advance/recede by a rate of 2.5 /orbit, 
allows for rendezvous in less than seven days. The ascent 
stage is not thermally sensitive at this point in the 
mission, therefore the rendezvous passage presents no 
thermal problems. 


The final ascent trajectory has the following 

characteristics : 


Ascent Thrust (82 seconds) 4020.3 m/sec 

Delta V Launch in a / s 

Maximum Acceleration at Launch 61 > q* 

Launch Angle 

Max. Mach (75 km alt.) 2732. N/m2 

True Anomal^Traversed^ exoatmospheric/atraospheric ) 

Time in Atmosphere (up to 250 km altitude) 


Flight Path Angle at Insertion 

V (600.9 km Altitude) 

V Circular (MOV) 

Delta V insertion . . . . 

Delta V (combined launch & insertion) 

Maximum Acceleration at Insertion 


196 sec 
2.92° 

2456.1 m/s ec 
3271.4 m/sec 
828.0 m/sec 

4848.3 m/sec 
< 10 g's 


Total Delta V: (m/s) 

Positioning 
Deorbit 
Retrof ire 
Ascent 
Insertion 
Rendezvous 


8.1 
167.9 
42.9 
4020.3 
828.0 
8 . 1 


5075.3 
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Mission Overview Diagram 
Mission Profile 
Ideal Descent 
Equations 
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Detailed Analysis 

Force Balance Diagram 

Pertinent Equations 

iterative Time Stepped Procedure 
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Linear Plot 
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Graph 
Table 

Final Descent Profile (e = 0.10) 

Graphs 

Descent Altitude vs. 

Descent Velocity vs. 

Descent Velocity Vs. 

Ascent Analysis 

Minimum Insertion Delta v 
At Optimum Launch Angle 
Graph 

Chart , . 

Final Ascent Profile Spreadsheet Thrust 

Launch Delta V Due To Constant 1500 N Thrust 

Mach Number a vs! n Altitude for Ascent and Descent 

Chart (Includes Dynamic Pressure) 
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Delta V 
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3 
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4 
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5 
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162.115 

6 
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190.021 

7 

53.244 
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8 
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9 
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Section E: 

Specialized Structures 


Tim Gibson 



Parachute 


During the descent through the martian atmosphere, the lander 
will employ a parachute to deaccelerate to a velocity at which the 
liquid engine can prepare the lander for touchdown. The size of the 
parachute, measured by its reference area, the area of a circle 
whose diameter is equal to the distance across the fully deployed 
parachute, must be large enough to create enough drag so that the 
required final or "terminal" velocity will be attained, without too 
much deacceleration during the opening "jerk". It can be estimated 
from Air Force research (ref. 3) that the peak deacceleration during 
opening jerk is slightly more than 3 times that of the deacceleration 
calculated from a force-balance relationship. The lander is limited 
to a maximum acceleration of 10 earth g's. This limits maximum 
force-balance calculated deacceleration to under 3.33 g's. 

The force-balance relationship is: 





Oinam-i. i 

OF PQUri * V ;oA.U'7 



which yields: D - mg = ma 


where: drag, D = C D V 2 3 v 2 S 

drag coefficient, C D 
atmosphere density, 3 
velocity, V 
reference area, S 
mass of the lander, m 
deacceleration, a 

so, we have: C D V 2 3 v 2 S - mg = ma 

From this, by setting a = 0, we can obtain the reference area 
required for a desired terminal velocity. 

S = 2 m 9 / c D 3 V 2 

or, the maximum force-balance deacceleration. 

a = ( CD 1/2 3 v 2 S - mg )/ m 

where V is the deployment velocity. 

By starting at V = deployment velocity, and recalculating V 
every time interval At by numerically integrating the acceleration 
(see Appendix j£A ). the maximum descent distance and time to reach 
terminal velocity can be determined. This method yields a maximum 
because the actual acceleration rate is, as mentioned earlier, higher 
than the force-balance method predicts. However, the maximum is 
desirable for determining the minimum altitude at which the 
parachute can be deployed safely. Low altitude parachute 
deployment is desirable to minimize the effect of cross winds. 

Initially, the lander was to use a 2-stage parachute. The small 
first stage, 3.33 m diameter, limited peak deacceleration to 3.14 gs 
at a deployment velocity of 280.00 m/sec at 1 km altitude (fig. /j. , )• 

The larger second stage, 9.9 m in diameter, limited peak 
deacceleration to 2.69 g's at a deployment velocity of 88.00 m/sec 
(fig £ £ ) and attained a terminal velocity of less than 30.00 m/sec 


in less than 75 m. This configuration could be deployed at less than 
1 km above ground level and maintain a 3 times maximum g's value 

under our lOg mission limit. _ 

The parachute is an extended skirt canopy type (fig J=. 3 ) which 
has a higher drag than flat circular or cross-form types but has a 
longer opening time and correspondingly has lower peak opening 
forces The canopy is fabricated of 2.60 ounce per square yard 
(89 9/ m 2 ) dacron cloth, per MIL-C-7350, type 1. The suspension 

lines are tubular braided dacron with an ultimate tensile strength of 
2740 N. The approximate weight of the canopy can be calculated 
from the surface area of 1/2 of a sphere: 

^surface = 2 K r 2 

and a required storage volume (chute pack) based on a maximum 
packing density of = 600 ^9/ m 3 . 

So, the mass of the chute canopy, m canopy = 89 A sur(ace = 89 (2 it r 2 ), 

and the volume of the chute pack, v pack = ( 89 ) ( 2 n r ^ / 600x1 0 3 - 

This mass and volume dependence on the radius squared led to 
an optimization where the weight and size of the parachute were 
compared to the weight of added liquid fuel to accomplish the same 
deacceleration. It was determined that decreasing the parachute 
from 9.9 meters in diameter ( m canopy = 13.7 kg ) to 4 meters in 
diameter ( m canopy = 0.75 kg ), while adding = 4 kg of fuel, would 
decrease mass by - 8.9 kg. Also, the volume required for the chute 

pack went from 0.0228 m 3 to 0.00125 m 3 . 

Using the 4 m diameter parachute, an analysis was done 
(fig. E ) for the earliest possible deployment, 5 km above ground 
level, with the highest deployment velocity possible, V = 248.54 
m/sec Maximum force-balance deacceleration was 2.47 g's, with a 
distance to terminal velocity of - 0.5 km. Next, an analysis was 
done (fig. ^ S') for a late deployment. -3 km altitude, with a 
corresponding terminal velocity of 187.4 m/sec. Maximum force- 


balance deacceleration was 2.48 g’s, with a distance to terminal 
velocity of = 325 m. 
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c 


c 

50 


c 


60 


pr ogra m mars 
implicit none 
real*4 cd, rho, s, 
real *4 del tat 
integer^ i 


g t m , y, y term , t , h , 


cd = 1.0 
del tat = 0. 1 
r+io = 0.0127 
s = 12.5664 
g = 3.73 
■ = 100.0 
v = 187.4 

vterm = sqrt( (2 . 0*m*g )/(cd'* , rtio’*'s ) ) 


»rite(6,50) uterm 

format< lx, ’v term = f7.3, m/sec ,/, lx, 

* 5x,h (meters)' ,5x, accel <g s) , 

* 5x, 1 time (sec)’ ) 


t = 0.0 
h = 0.0 

do 100 i = 1, 200 
t = t + del tat 

a = <0 . 5 * c d * rho * <v** 2 )*s - m*g)/m 
h = h + v+del tat - 0.5*a*(del tat**2> 
onem = sqrt((v**2> - 2. 0 * 0 * 6 ) 


accel = a/9.81 
«ri te<6,60) v, 
format< lx, f7.3, 
if ( v . le. (1. 


h, accel, t 

7x, f7.2, 10x,f7.5,8x, f5. I 
0000 1*w term)) go to 200 


o ~ vnem 
enddo 
cont i nue 
end 


i, accel 


(m/sec)’ , 


200 
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V (m/sec) 

h (meters) 

280.000 

27.85 

276.911 

55.39 

270.816 

82.32 

261 .952 

108.39 

250.698 

133.34 

237.537 

156.98 

223.025 

179.19 

207.745 

199.89 

192.264 

219.05 

177.101 

236.71 

162.700 

252.94 

149.412 

267.85 

137.489 

281 .57 

127.079 

294.26 

118.233 

306.07 

110.917 

317.15 

105.024 

327.65 

100.398 

337.68 

96.855 

347.37 

94.202 

356.78 

92.258 

366.01 

90.861 

375.09 

89.874 

384.08 

89.189 

393.00 

88.720 

401 .87 

88.403 

410.71 

88.193 

419.53 

88.055 

428.33 

87.966 

437.13 

87.909 

445.92 

87 . 873 

454.71 

87.850 

463.49 

87.836 

472.28 

87.828 

481 .06 

87.823 

489.84 

87.820 

498.62 

87.818 

507.41 

87.817 

516.19 

87.817 

524.97 


cel (g s) 

time (s^ 

3.14891 

. 1 

3.07227 

.2 

2.92355 

.3 

2.71320 

.4 

2.45618 

.5 

2.16995 

.6 

1 .87222 

.7 

1 .57901 

.8 

1 .30315 

.9 

1 .05366 

1 .0 

.83567 

1.1 

.65093 

1 .2 

.49855 

1 .3 

.37586 

1 .4 

.27919 

1 .5 

.20451 

1 .6 

.14782 

1 .7 

.10549 

1 .8 

.07435 

1 .9 

.05178 

2.0 

.03563 

2.1 

.02423 

2.2 

.01629 

2.3 

.01082 

2.4 

.00710 

2.5 

.00461 

2.6 

.00295 

2.7 

.00187 

2.8 

.00117 

2.9 

.00072 

3.0 

.00044 

3.1 

.00027 

3.2 

.00016 

3.3 

.00009 

3.4 

.00005 

3.5 

.00003 

3.6 

.00002 

3.7 

.00001 

3.8 

.00001 

3.9 
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V (m/sec; 

h (meters) 

88.000 

8.67 

85.356 

17.08 

80.265 

25.00 

73.286 

32.24 

65.22? 

38.70 

56.989 

44.35 

49.372 

49.26 

42.953 

53.54 

38.017 

57.33 

34.550 

60.78 

32.322 

64.00 

31 .004 

67.10 

30.27? 

70.13 

29.907 

73.12 

29.727 

76.0? 

29.644 

79.06 

29.60? 

82.02 

29.594 

84.98 

29.58? 

87.94 

29.587 

90.90 

29.586 

93.85 


cel ( g s ) 

time (sec) 

2.69566 

. 1 

2.51574 

.2 

2.18488 

.3 

1.76431 

.4 

1 .32630 

.5 

.93109 

.6 

.61311 

.7 

.38056 

.8 

.22368 

.9 

.12494 

1 .0 

.06647 

1 .1 

.03371 

1.2 

.01629 

1 .3 

.00749 

1 .4 

.00328 

1 .5 

.00136 

1 .6 

.00053 

1.7 

.00020 

1 .8 

.00007 

1 .9 

.00002 

2.0 

.00001 

2.1 
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y term = 

90 .802m/sec 

V (m/sec> 

h (meter 

248.540 

24.73 

246.118 

49.23 

241 .337 

73.25 

234.372 

96.58 

225.501 

119.03 

215.086 

140.45 

203.546 

160.73 

191 .323 

179.80 

178.860 

197.63 

166.568 

214.25 

154.806 

229.69 

143.868 

244.05 

133.971 

257.43 

125.253 

269.93 

117.773 

281.70 

111.520 

292.84 

106.424 

303.48 

102.372 

313.71 

99.226 

323.63 

96.835 

333.31 

95.056 

342.81 

93.757 

352.19 

92.826 

361 .47 

92.169 

370.69 

91 .713 

379.86 

91 .401 

389.00 

91.191 

398.12 

91 .051 

407.22 

90.960 

416.32 

90.900 

425.41 

90.863 

434.49 

90.839 

443.58 

90.824 

452.66 

90.815 

461 .74 

90.810 

470.82 

90.807 

479.90 

90.805 

488.98 

90.804 

498.06 

90.803 

507.14 


cel ( g s > 

time (sec) 

2.46841 

.1 

2.41317 

.2 

2.30569 

.3 

2.15289 

.4 

1 .96476 

.5 

1 .75316 

.6 

1 .53037 

.7 

1.30781 

.8 

1 .09505 

.9 

.89924 

1 .0 

.72492 

1.1 

.57427 

1.2 

.44746 

1 .3 

.34324 

1 .4 

.25941 

1.5 

.19329 

1 .6 

.14208 

1 .7 

.10307 

1 .8 

.07381 

1 .9 

.05220 

2.0 

.03646 

2.1 

.02515 

2.2 

.01714 

2.3 

.01153 

2.4 

.00766 

2.5 

.00503 

2.6 

.00326 

2.7 

.00209 

2.8 

.00132 

2.9 

.00082 

3.0 

.00050 

3.1 

.00031 

3.2 

.00018 

3.3 

.00011 

3.4 

.00006 

3.5 

.00004 

3.6 

.00002 

3.7 

.00001 

3.8 

.00001 

3.9 
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\> term = 

68 . 369m/sec 

V (m/sec) 

h (meters) 

187.400 

18.62 

184.971 

37.00 

180.194 

54.91 

173.304 

72.14 

1 64 . 666 

88.51 

154.741 

103.91 

144.04<S 

118.25 

133.100 

131 .51 

122.387 

143.71 

112.318 

154.91 

103.209 

165.20 

95.269 

174.71 

88.593 

183.56 

83.180 

191 .87 

78.940 

199.76 

75.731 

207.32 

73.377 

214.66 

71 .701 

221 .83 

70.539 

228.88 

69.755 

235.85 

69.236 

242.78 

68.902 

249.67 

68.690 

256.54 

68.559 

263.39 

68.479 

270.24 

68.432 

277.08 

68.404 

283.92 

68.389 

290.76 

68.380 

297.60 

68.375 

304.44 

68.372 

311 .29 

68.371 

318.11 

68.370 

324.95 


accel <g s) 

time <se 

2.47641 

. 1 

2.40283 

.2 

2.26095 

.3 

2.06283 

.4 

1 .82535 

.5 

1 .56750 

.6 

1 .30756 

.7 

1 .06080 

.8 

.83816 

.9 

.64593 

1 .0 

.48624 

1 .1 

.35805 

1 .2 

.25821 

1 .3 

.18257 

1.4 

. 12667 

1 .5 

.08628 

1 .6 

.05774 

1.7 

.03796 

1 .8 

.02452 

1 .9 

.01556 

2.0 

.00970 

2.1 

.00594 

2.2 

.00358 

2.3 

.00211 

2.4 

.00122 

2.5 

.00070 

2.6 

.00039 

2.7 

.00021 

2.8 

.00011 

2.9 

.00006 

3.0 

.00003 

3 . 1 

.00002 

3.2 

.00001 

3.3 



Heat Shield Structure 


The structure of the heat shield serves to hold the insulation 
material in position, and serve as an aero-braking device. The 1.5 
meter diameter heat shield is designed to withstand a maximum 
total aero-force of 6.9 kN which occurs at an altitude of 31 km and a 
corresponding free-stream velocity of 1985 m/sec. The structure in 
contact with the insulation encounters a maximum temperature of 

700 F for approximately 30 seconds. 

Due to its high stiffness-to-weight ratio and resistance to 
thermal effects, a beryllium honeycomb bonded with high- 
temperature thermoplastic adhesive was chosen for the heat shield 
structure. Four beryllium tubular braces serve to brace the 
structure (see fig. E (o )• 
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In order to conservatively size the structural components 
simplified models (see fig. E7 ) were chosen which would provide an 
adequate safety margin and facilitate analysis. 
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where: unsupported length, L 

load per unit length, P 
depth of honeycomb, h 
thickness of face sheet, t 
width of beam, w 

From strength of materials the stress is: stress = My /| 

now, y = 1 / 2 d + t , but since t« 1 / 2 h, y= 1 / 2 d 

and, by the parallel axis theorem I = 2 / 2 wt 3 + 2wt( 1 / 2 d) 2 
since 2 / 3 wt 3 « 2wt( 1 / 2 d) 2 , we have I = 2wt(V 2 d) 2 = 1 / 2 wtd 2 


-V 7 • -Jj j V 



- M( ' /2)h '('/ 2 )Wth2- M/ Wtd 


so, stress 


or, letting M = moment per unit width, stress = ^/td 

Also from strength of materials we know the maximum moment, M, 
of our model beams (fig. /r*7) f° r case 1 is equal to 

M (max.) = P*- 2 /8 . 

and for case 2 

M (max) = ^ 2 /8 • 

Using a maximum unsupported length of L = 1 m for case 1 and L 
0.5 m for case 2, we have: 

M (max) = Wq 

Now stress crit for beryllium = stress yjeld = 3.24x10® N / m and the 
distributed load, 

p = Aero-force/ area = (6.9 Kn/^ 75m) 2)(1m) = 3905 N / m 

Rearranging, we have: t d = P ^8(stress y ) = 3905/ 8(3. 244x1 0 8 ) 

=1.507x10-® 

Taking the minimum thickness for the beryllium face sheets to 

be t = 0.010" = 0.000254 m, => d = 1 507x1 0 6 /. 000254 

= 0.00593m = 0.23 in 


so, we will specify a beryllium honeycomb with a core depth, 
d = 0.25" = 0.00635 m. Checking, we have 


stress crit = 3905 /8(.00635)(. 000254) = 3.03x10® N / m 2 


which is over 6% under stress y 
Now, in order to determine 
density of the honeycomb core is 
estimating honeycomb 


the mass of the structure the 
needed. A "rule of thumb" for 


core density is 


^core = 2 ^ Be ^ max/, stress 


Be shear 


where: density of Be, 3 Be = 1826 k 9/ m 3 

shear stress of Be, stress 5 = 0.6 stress y 
and the maximum shear stress, again from strength of materials 


V m ,x -V- 1-982-5 N 

so, p„ re = 2(1 ®2 6 )(1 9 52 - 5 )/q 6(3 24 x 1 0 8 ) - 0 0376 

The lightest core material available is approximately 1.5 Ibs/ft or 
24 kg/m. Apparently the shear force is too low in this application 
for the "rule of thumb" approximation to be accurate. It does 
indicate though that it is safe to use the lightweight core. So, the 
mass of the structure is: 

Mass = 23 Be (7tr 2 t) + 3 core (*r 2 d) 


= nr 2 (23 Be t + 3 core d) 

= 7 i(. 75 ) 2 [ 2 ( 1826 )(. 000254) + 24(. 00635)] = 1.91 kg 

The model for the tubular braces is simple compression 
(see fig £ 8 )• The force in each tube may be as high as 


F = (1.5)( max - aero-force /# Qf tubes ) 



where: max. aero-force = 6.9 kN 

# of tubes = 4 

so F = ( 15 ) ( 690 °)/ 4 = 2588 N 



where: diameter, d 

radius, r 
length, L 
wall thickness, t 
force, F 

For a thin wall tube under compression there are three design 
constraints: 

1) stress limit: stress = ^dt 

2) Euler buckling: F = k2 ^/\_ 2 


3) crippling: stress = 0.3EV r 



For a first iteration, the stress limit and Euler buckling are solved 
together in order to find the optimum tube dimensions, where 

I = m^t for a thin-wall tube. 

Now, solving 1) and 2) for the force, F, yields 

(stress)27trt = Tt2 ^ 7rl ' 3 Vi _2 

so r- 2 (stress)L 2 /^ 2 E = 2(.652) 2 (3.24x1 0 8 )/ tc2(2 8g6x1 Ql 1} 

= 0.009817 m. 

This is a tube with a diameter of 0.0196 m or 0.773". Solving for 
the wall thickness, t = 0.0001313 m = 0.0052 in. In order to specify 
a tube which might be commonly available, we will use 0.75 in. 
(0.01875 m) diameter and 0.010 in. (0.000254 m) wall thickness 
beryllium tubing. Checking, we have: 

1) stress limit, F = (stress)2nrt = 4848 N 

2) Euler buckling, F = re2 ^/|_2 = 4421 N 

3) crippling, stress = 0.3EV r = 23.5x10 8 

which are all considerably over any anticipated forces. Now, the 
mass of 

each tube is, mass/tube = 7tdtL3 Be 

= rc(. 01 875) (.000254) (.652) (1 826) = .0178 kg 
So, the total for four tubes is 4(0.0178) = .0712 kg 


Therefore, the mass of the entire structure is: 1.91 + 0712 

kg 


1.98 


Soil Sample Collection 


The mission requirement is to bring 1 kg of Martian soil back 
to the orbiting MOV (Mars Orbiting Vehicle). This will be 
accomplished with the use of an extendable surface sampler, which 
consists of a collector head attached to the end of a retractable 
boom. The boom is constructed from two ribbons of beryllium 
welded together along the edges. When extended, the two layers 
opened to form a rigid tube. When retracted, the boom flattens onto 
a storage spool. A flat cable sandwiched between the boom layers 
transmits electrical power to the collector head. 

The collector head is a scoop with a movable lid. To fill the 
scoop, the boom is extended along or into the surface. Once full, the 
lid closes. The boom is then elevated such that the scoop can be 
dumped into a container at the top of the lander (fig .£?)■ The nose 
cap at the top of the lander's aero-shell (fig. £ /O ) is moved upward 
by a threaded driveshaft and, when free of the aero-shell, is then 
rotated out the way for dumping into the sample container. 

Design considerations include tube buckling and bending 
stresses, and torque for pushing the scoop into the Martian soil, 
lifting the sample into the container, and opening and closing the 

container. 

Data . Martian soil: 

density, 3 = 0.7 to 3.2 9/ cm 3 . 

N 2 

maximum penetration resistance, PR = 6 / cm C m 
Scoop 

Volume required for 1 kg sample, 3 = 0.7 9/ cm 3 
1 / 7 00 = 0.00143 m 3 = 1430 cm 3 



Jfc co P 



where: s = scoop side, L = scoop length, t = material thickness. 

Since the collector is being driven into the soil, titanium is used as 
a material for the scoop due to its toughness. The worm-gear motor 
should have enough torque to close the lid against maximum 
penetration resistance. Now, the force needed to drive the scoop 
into the soil, 
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F = PR (length) (scoop cross sect, area) 

where: scoop area = 4 st, and length, L = 16.5 cm, s = 10.16 cm, 
and t = 0.010 in = 0.0254 cm, 

so, F = (PR) (L) [4 (s) (t)] = 6 (16.51) [4 (10.16) (0.0254)] 

= 102.26 N 

A force of F = 200 N will be used in the design of the sampler 
apparatus to insure an adequate safety margin. 

Mas! 

The dimensions of the mast (boom) will be governed by 
maximum compressive stress, Euler buckling, and bending. The 
bending requirement will be to lift the scoop, sample, and mast into 
the "dump position." 

Maximum length of mast = 1.5m 



The wall thickness of the beryllium mast is taken to be 0.0508mm, 
which is common for the welded seam type masts. 

t = 0.0508 mm = 5.08 x 10‘ 5 m 
= 0.002 " 

In order to find the minimum diameter, maximum compressive 
stress is considered: 


stresSy = F //\ => F = stress yje | d A - stress y 7idt 

=> d = F / stre ss y nt = 200/ (3.24x10 8 ) it (5.08x10' 5 } 
- 0.00387 m 

SO, d = 0.387 cm 

7l 2 EI/, p _7i 2 E7cr 3 t/. 2 

Now consider Euler Buckling: F = L L 

so, r = ( FL2 / 7t 3Et) 1/3 = (200(1 - 5 ) 2 / n 3 (2. 896x1 0 M )(5. 08x1 O' 5 )) 

= 0.00995 m 


d = 0.01 991 m = 1 .991 cm = 0.784 " 
and finally bending, from strength of materials: stress = ^7/| , 
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Mass of Be mast of length 1.5m, 

rn mast = rrdLta Be = jcd(1 ,5)(5.08x10- 5 )(1826) = 0.437d 

Now, the moment, M = m mast ( 1 / 2 L) + m samp | e L + m scoop L. 


IS 







worm gear motor 



so M = ,437d(.75) + (1.5) + .5288(1.5) = 3277d + 2.2932 
and I = nrH = { k / q ) 6 3 \ 

so stress y = M V/| = (- 3277d + 2 - 2932)4 /jid 2 (5.08x1 O' 5 ) 

= (1'31d + 9 - 1728 )/(i .596x10' 4 d 2 ) 
but stress y = 3.24x10 s , so, rearranging, we have 
5.171x10 4 d 2 - 1 ,31d - 9.1728 = 0 


or d 2 - 2.533x10' 5 d - 

1.774x1 0‘ 4 = 0 

solving for d: 

d 

F(dl 

.01991 

.0002185 

.05 

.00232 

.01 

-.000077653 

.015 

.00004722 


so, for bending 0.010m < d < .015m. The Euler buckling in 
compression is the limiting case, therefore the diameter is chosen 
to be d = .01991 m. 


Finally, an analysis of crippling failure: 



Fcri p/ TC dt = °- 3E = 2t/ d 

F crip = 0.3rc(2.896x10")2(5.08x10' 5 ) 2 = 1408 N 

The apparatus is designed for a maximum of 200N, so there is no 
problem here. 

Reel Assembly 




Worst case torque required would occur while forcing scoop into 
soil, so T = Fr where F = 200 N 

r = 7.25 cm = 0.0725 m 
T = 200(.0725) = 14.5 Nm 
= 128.3 in lbs 

A Hurst Instrument Motor #2602-001 has a mass of 0.68 kg and 
delivers 250 in lbs of torque at 1 RPM, which corresponds to a scoop 
speed of 1 .5 cm / sec The maximum power requirement of each motor, 

P = 54 watts, so the electrical power requirements per sample 
collection are: 

time for unreeling at 1.5 cm /sec ~ ^0 sec, 2 (100) = 200 sec, 

motors running = 3, => 32400 watt-sec, 

time to rotate through boom swing = 30 sec, 2 (30) = 60sec, 

motors running = 3, => 9720 watt-sec, 

time to rotate to dump position at 1 RPM = 15 sec, 2 (15) = 30 sec, 

motors running = 4, => 6480 watt-sec, 

time to dump * 15 sec, 

motors running = 4, => 3240 watt-sec. 

So, the total energy required per try = 51840 watt-sec. 

The lander will have 300 watts available for 10 minutes (600 sec) = 
180000 watt-sec. Therefore there will easily be enough energy for 
at least 3 complete attempts. 
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CONCLUSIONS 


We believe our conceptual design to represent the best version 
of the three being presented this quarter. This design has by far the 
lightest weight and shows much evidence of the optimization 
proceedures we used to improve our designs while the competing 
designs were busy attempting to incorperate some of the best 
features of our earlier models. Although our design is by no means 
perfect and optimized in every possible way, it certainly represents 
a very solid foundation upon which a final design may eventually be 
built. Our team has enjoyed this project very much and we all feel 
we have learned a great deal while working on it and we look 
forward to working on similiar challenging projects in the world 

outside of UCLA. 


